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TWO- AND THREE-DIMENSIONAL UNSTEADY LIFT PROBLEMS 

IN HIGH-SPEED FLIGHT 1 

By Harvard Lomax, Max. A. Heaslet, Frankly** B. Fuller, and Loma Sluder 


SUMMARY 

The problem, of transient lift on two- and three-dimensional 
wings flying at high, speeds is discussed as a boundary-calue 
problem for the classical leave equation. Kirchhoffs formula 
is applied so that the analysis is reduced, just as in the steady 
state, to an investigation of sources and doublets. The appli- 
cations include the evaluation of indicial lift and pitching- 
moment curves for two-dimensional sinking and pitching 
wings flying at Mach numbers equal to 0, OH, 1.0, 1£, and 
2.0. Results for the sinking case are also given for a Mach 
number of 0.5. In addition, the indicial functions for super- 
stmic-edged triang^ilar wings in both forward and reversed flow 
are presented and compared with the two-dimensional values. 

INTRODUCTION 

The usual idealizations introduced in the development of 
linearized aerodynamic theory describe a frictionless, per- 
fectly elastic, model fluid. As is well known, the effect 
of small disturbances in such a fluid can be analyzed by 
means of the familiar wave equation. 

The solution to the wave equation known as Kirchhoff's 
formula (see reference 1) is found to be of considerable use in 
unsteady-motion problems involving thin wings with super- 
sonic edges. The problem is reduced to one of summing 
elementary solutions, analogous to sources and doublets 
in steady flow, over a region determined by the position of the 
wing as well as its traversed path. The theoretical develop- 
ment leads naturally to the concepts (defined later) of, first, 
inverse sound waves, which have a counterpart in the Mach 
forecones used in steady-state wing theory; second, acoustic 
plan forms; and, third, homogeneous flow, which reduces 
in part to the familiar conical flow as the wing approaches 
a steady supersonic velocity. 

There are several simple types of unsteady motion on 
which the analysis can be based. The so-called indicial 
motion, in which the velocity undergoes an initial discon- 
tinuous change will be considered here. It is possible to 
conceive the perturbation field due to the unsteady motion 
in two slightly different ways. For one, it can be supposed 
that the wing has been traveling at the constant velocity 
r 0 for an in finitely long time and then, at time equals zero, 
starts suddenly to sink without pitching (or pitch at a 


constant angle of attack) while maintaining the forward 
velocity TV On the other hand, the wing may be considered 
to be at rest in still air until at time zero it starts suddenly 
either to sink or to pitch and, at the same instant, attains 
the forward velocity T" 0 . The latter physical picture will be 
used in this report. Problems of unsteady motion can also 
be approached with the initial assumption that the velocity 
potential depends harmonically on the time. These two 
approaches are quite compatible in that they can be related 
through the use of superposition methods (Duhamel's 
integral, Fourier’s integral) of the operational calculus. 

Detailed results for two classes of indicial responses (time 
response to a step input) will be given in this report. The 
first of these is the complete set of responses {c ta , c„ a , c tf , c„ f ) 
for two-dimensional wings flying at Mach numbers equal 
to 0, 0.8, 1.0, 1.2, and 2.0. In addition, results for c, a and 
c ma are given for a Mach number equal to 0.5. The part of 
the analysis pertaining to the response of wings traveling at 
subsonic speeds is lengthy and somewhat tedious regardless 
of the method of approach. With the use of indicial func- 
tions, however, the calculations are reasonably straightfor- 
ward, especially for Mach numbers around O.ff to 1.0. 
Further, the use of indicial functions sheds considerable 
light on the manner in which Mach number variations 
affect the section aerodynamic characteristics. The results 
for the two- dim ensional wings traveling at supersonic speeds 
are compared with the indicial responses developed by 
rectangular wings of aspect ratios 2, 4, and 6, the latter 
curves having been' presented by J. W. Miles in reference 2. 

The second class of indicial responses considered is that 
concerned with the forces and moments induced on sinking 
and pitching triangular wings. First, the loading on a flat 
triangular wing with supersonic edges undergoing an indicial 
sinking motion is determined. (See also reference 3.) Then 
a simplified method is developed whereby total lift and 
pitching-moment coefficients for the wing with supersonic 
edges may be obtained. (See also reference 4.) Again a 
complete set of indicial responses is presented for supersonic- 
edged wings in both forward and reversed flow. Lastly, 
the triangular wing with subsonic edges is partially analyzed, 
and the indicial responses for a slender triangular wing are 
given. 


i Supersedes N'AC A TM 2103, “The Indicial Lift and Pitching Moment for a Sinking or Pitching Two-Dimensional Wing Flying at Subsonic or Supersonic Speeds'’ by Harvard Lomax 
Mai, A. Heaslet, and Loma Slnder, 1951; MAC A. TX 23S7, “Three-Dimensional Unsteady Lift Problems In High-Speed Flight— The Triangular Wing” by Harvard Lomai, Mai- A, Heaslet* 
and Franklyn B. Fuller, 1951; and also contains material from NACA TX 2256, ‘‘Three-Dimensional, Unsteady-Lift Problems In High-Speed Flight— Basic Concepts** by Harvard Lomax, 
Max. A. Heaslet, and Frr.nklyn B. Fuller, L95L 
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SYMBOLS 

speed of sound in the undisturbed fluid 
chord of wing 

lift coefficient — \ 


/. J*L_\ 

\l PoVaS J 


indicial lift coefficient due to angle-of-attack 

change, without pitching ^ i a == ^~| ^ 

indicial lift coefficient due to pitching for a 
whig rotating about its leading edge or apex 




pitching-moment coefficient, positive when trail- 
ing edge tends to sink relative to leading edge 

( mo ment \ 

~ P«TVc5 ) 

indicial pitching-moment coefficient due to angle- 
of-attack change (without pitching) measured 
about the leading edge or apex; positive when 
trailing edge is forced downward with respect 

to the leading edge or apex ) 

indicial pitching-moment coefficient due to pitch- 
ing measured about the leading edge or apex, 
for a wing rotating about its leading edge or 
apex; positive when trailing edgo is forced 
downward with respect to the leading edge or 

bc m 

L KfoX-oJ 

two-dimensional lift coefficient / y- \ 

\2 poF °* c / 


two-dimensional pitching-moment coefficient 


( moment \ 

3 '• w) 


cotangent of sweep angle (cot A) 

Mach number in the undisturbed fluid 

M 0 f Ap , 


Ml f iS-dy 

2as ,/-« Qo 
rnMgVo C‘ A P j 


loading coefficient (pressure on the lower surface 
minus pressure on the upper surface divided by 
free-stream dynamic pressure) 

free-stream dynamic pressure (~ p 0 F 0 2 ^ 

dimensionless rate of pitclmig 

local semispan of wing 

whig area 

time 

04\ Gto t'fc 


u, v, mo perturbation velocities hi x, y , z directions, re- 
spectively 

V 0 free-stream velocity 

x, y, z Cartesian coordinates, fixed relative to the fluid 
at infin ity 

(a:/c)e. 5 . distance of center of pressure from leading edge 
or apex, in percent chord 

a angle of attack (angle between flight path and 

plane of wing), radians 
da 

“ IF 

$ 

r circulation 

5(f') Dirac S function, normalized with respect to V, 

thus 5(t')=0 for t'j& 0, 5(0) = °°, and 

f S(f')dt ’= 1 

— oo 

A discontinuity in the quantity hi question across 

the plane of a wing 

6 wing angle of pitch relative to horizontal, posi- 

tive when trailing edge lies below leading edge, 
radians 

6 wing rate of pitch 

A angle of sweep of leading edge, positive for sweep- 

back 

Po free-stream density 

r 0 ■ chord lengths traveled 

<p perturbation velocity potential 

SUBSCRIPTS 

u upper surface of a wing 

/ lower surface of a wing 

x, y, z, t indicate differentiation with respect to the vari- 
able in question 

PART I — THE USE OF INDICIAL FUNCTIONS 
IN UNSTEADY LIFT PROBLEMS 

In the first part of this report, a discussion of three un- 
steady lift problems will be given. The three problems are: 

(i) Determination of the indicial response in lift and 
moment on a sinking or rotating wing; 

(ii) Determination of frequency response in lift and mo- 
ment on a fluttering wing; and 

(iii) Determination of frequency response in lift and mo- 
ment on a slowly oscillating wing. 

The method by which the latter two problems may be 
solved with the aid of the solution to the first is also de- 
scribed, as well as the application of these results to the de- 
termination of lift and moment on wings undergoing arbitrary 
maneuvers. 

THE INDICIAL FUNCTIONS 

By definition, an indicial function is the response to a 
disturbance which is applied abruptly at time zero and is 
held constant thereafter; that is, a disturbance given by a 
step function. For example, if the angle of attack of a wing 
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varies with time as shown in figure 1, the corresponding re- 
sponse, also shown in figure 1, is designated as the indicial 
lift coefficient due to angle of attack. Four such indicial func- 
tions will be evaluated, namely, Cl*, Cm/, Cl/, Cm/. The 
primes on the coefficients indicate that the axes about which 
pitching motion occurs and pitching moments are measured 
either coincide with the leading edge of the wing or with a 
line through the apex normal to the root chord of the wing. 

The equations which transform these functions to those 
for a wing pitching about an axis a distance ac back from the 
leading edge and having its moment center a distance be back 
from the leading edge are simply 


(, L a =C La 

<\ a =C m /+b C La 
r L =CL/-a C La 

r-C' + b C L '—a C'-ab Cr 


( 1 ) 




In linearized, thin-airfoil theory the boundary condition 
that applies to the indicial functions due to angle of attack 
is simply 

w a =w t = — Voa (2) 

over a certain planar area in the xyz space at the time t'. 
If we consider a coordinate system fixed relative to the fluid 
at infini ty and a wing with span measured along the y axis, 
moving away from the origin along the negative x axis with 
a velocity T’ 0 , then equation (21 applies to the area in the 
g=0 plane occupied by the wing at any given time. It is 
further required that the <? be continuous everywhere except 
across the wing plan form and wake. In the case of the 
indicial functions due to pitch for a wing rotating about its 


leading edge or apex, the boundary conditions for a wing the 
leading edge or apex of which is at the origin at f'= 0 are 
that the upwash be given by the expression 


w u =Wi= — (x+Vot') 6 (3) 

over the same region in the 2=0 plane as for the angle-of- 
attack case and, again, that- <p be continuous except across 
the wing and its wake. The angle of pitch, 9, is taken as 
positive when the trailing edge is lower than the leading edge, 
and 6 is the time derivative, dd/dt', positive when the trailing 
edge is falling with reference to the leading edge. 

The difference between 8 and a is illustrated in figure 2 (a). 
The angle of attack a is the angle between the flat wing 
surface and the tangent to the flight path of some point 
fixed at a distance ac back from the wing leading edge. For 
example, in applications involving the dynamic behavior of 
an entire airplane the distance ac would usually be taken as 
the distance back to the center of gravity of the airplane. 
The angle 8 is the angle between the flat wing surface and 
the horizontal. Figure 2 (b) shows a wing undergoing a 
sinusoidal angle-of-attack variation with a zero angle of 


Direction of wing motion 




(b; 



(a) Angle convention. 

Cb) Sinking wing. 

(c) Pitching wing. 

Fig chi 1— Definition, of angles of attack and pitch. 
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pitch throughout. For convenience and in order to dis- 
tinguish from the pitching wing, such a wing will be referred 
to as a sinking wing. Figure 2 (c) shows a wing undergoing 
a sinusoidal angle-of-pitch variation taken about the axis ac 
back from the leading edge. In the case shown the wing has 
a constant (zero) angle of attack. 

The variable q to which the lift and pitching-moment 
coefficients of the pitching wing are referred is equal to 
QcjV, o, a dimensionless form of expression for rate of pitch. 


£ -+ i£ -"S : ij jfw£ 

R.'+Uf.'-ta 

lim Vp-arf K«'-«+ 


ADAPTATION OF THE INDICIAL FUNCTIONS TO MANEUVERS 

Since the theory is linear, the lift and moment on a wing 
undergoing an arbitrary maneuver (in which both a and 8 
remain small) can be calculated from the indicia! functions 
by the principle of superposition. For example, consider a 
wing performing some arbitrary maneuver involving small 
variations of both a and 6. Let a and 9 be defined * with 
respect to the flight path of a point ac back from the leading 
edge. The appropriate indicial functions can then be cal- 
culated from those referred to the flight path of the leading 
edge by means of equations (1). (The moment center is 
still arbitrary at a point a distance be back from the leading 
edge. However, for convenience, b is usually taken to be 
equal to a.) Finally, by Duhamel’s integral, the lift and 
moment induced on the wing by the arbitrary maneuver are 
given by the equations 


f‘\c La (f'-V) « ft') + C Ls (f'-hO q ft')] du 
Cm= Mo lCm ° (*' -*<')« ft') + ft')l dW 


Y (4a) 


W+iW-Bm &£\c..’V-V>+ 

where u is the frequency and k is the reduced frequency, 
oic/2V. The primes on the quantities again indicate that 
the wing is pitching about and the moments are measured 
about the leading edge of the wing. 

The responses in lift and moment for harmonically oscil- 
lating wings can be expressed in terms of their absolute 
values and phase shifts with reference to the period of the 
forcing disturbance. Thus, if an angle of attack variation 
given by the equation a=c i “ t ’ is impressed on the wing, the 
response in lift would be [c + In terms of the 
indicial functions this can be written 

| Ci f e, 

a dt Jo a 

which, by comparison with equations (5), can be re-expressed 
in terms of the flutter derivatives as 


or alternatively 


c .Je<*=4 k(I, t -iLd = 4 it V£ s 2 +£»*e^ ,an 




y (4b) 


where, for example, Cijjf—ti) means that the indicial 
function (? L is to be evaluated at the time 

«Jf 


THE FLUTTER DERIVATIVES 

By definition, a flutter derivative is the response to a 
harmonic oscillation in vertical displacement (angle'of attack) 
or pitch, such oscillation having continued so long that all 
transient effects of its origin are damped out and the induced 
forces and moments are periodic. In the notation used in 
reference 5 the induced responses are given in terms of the 
flutter derivatives L u Lz, L s ', L/ and Mi, M t ', M/, M t '. 
The relation between the indicial responses and the flutter 
derivatives follows from equations (4) and can be written 


Similarly, each of the stability derivatives can be expressed 
in terms of the indicial functions and the flutter derivatives, 


hence: 



Term 

Maximum value 

Phase angle 


4A-VT7+X? 

I 

lit 


2 kjMS'+Mt” 

r Ml 

-tan 

>■ 

i£(c ti +c,/) 

2£y 'W’+LS* 

lan "it 



- , Mi 

lan J!7 J 


where it is assumed that unit absolute values of angle of 
attack and pitch are impressed and the term a=c 

t a 


1 Notice that if o' and V ora the angle ot attack and pitch measured with respect to the flight path of the leading edge, and a and e the same angles measured with respect to ths flight path 
of a point ac back from the leading edge, the relation between the two sets of angles is given, for s mall deflections, by 
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A WING UNDERGOING SLOW HARMONIC OSCILLATIONS 

The boundary conditions for these problems are exactly 
the same as those used for evaluating the flutter derivatives, 
being in fact a special case of the latter when the frequency 
of oscillation is very low. This added simplification is im- 
portant, however, because the frequency of angle of attack 
and pitch oscillations for the entire airplane is, relative to 
the flutter frequencies, very low. Further, special methods 
based on the assumption that the wing is oscillating slowly 
have been devised. Thus, if the lift and moment are ex- 
panded in powers of the frequency, solutions for the coeffi- 
cients of the lowest-order terms have been found for both 
triangular and rectangular wings- flying at supersonic speeds. 
(See, e. g., references 6, 7, S, and 9.) 

For a given set of indicial functions these expansions can 
lie carried out with relative ease. Consider the case of L t 
and Lt as given by equations (5). If the indicial curve for 
r t ({') contains no pulse function, 3 the expression for Li+iLi 
can be written 


Li + iii=lim 

t ' — rOO 


ie 


-tvl* 


7( 0 £ 'M c °)- a 


where Ac f it') is the difference between the indicial lift and 
its asymptotic value, thus 

A<■< <I J') = <r< a (“)-C ^ „(^ , ) 


Since is zero, it follows 

Li + iLi= [c fa ( « ) - ia A c la {h')e -*•*»' cft/J 

or, in terms of the reduced frequency' parameter k=ucf2V 0 
and the number of chord lengths traveled T 0 =V 0 t'lc, 

L l J hiL i =^j' [c ta (~)-2ikj Q Ac j a (r 0 )e _2tl:ro </ r oj 


In case the flow is supersonic, these expressions can be 
evpanded in powers of k\ thus 

4^=2 | Act a (vo)dro — 4A 3 ro 3 ACj a (T 0 )d7-o+ - - - (7a) 

4l*i.=C( (<=) — 4Ar* f roAC( (ro)cfro-i- . - . (7b) 

,/o 


Similarly, the expressions for the other coefficients become 
— 2.l// = 2 f Ac^^Toldro— 4A ,S f ro 1 Ac s , o , (ro)cfro-f- . . . 

,/u */G 

(7 c) 

-2k I/ 1 '=c,;(®)-4i- ! I" r 0 Ac M J(r a )dr a + . . . (7d) 
“ Jo 

2A*L/=^ C{ o (eo) + [ACi^fro)— roACj a (ro)]di- 0 -l- . . . 

(7 e) 


2kL± =Ci'{a>) — J* Ae {a (ro)cfro— 2^*^ [2T 0 ACj t , (T a ) 

— ro s Ac, a (7-a)]dro+ ... (7f) 

-AW=| C K£t / (o3) + 2A- i J o ”[Ac„;(ro)-r 0 Ac s , tt '(r 0 )]dr 0 +. - - 

(7g> 

— r 0 *Ac» a '(r(i)]o?ro-|- - - - (7h) 


It is interesting to examine these equations briefly with 
regard to the problem of one-degree-of-freedom oscillatory 
instability. As was shown in reference 5, an oscillatory 
stability boundary is given by' the equation d/ 4 =0. Hence, 
a wing pivoted about its leading edge can be neutrally' stable 
if d//= 0. Such a condition arises only if the frequency is 
very low 4 for which case the above expansion for Ad// applies. 
Hence, to a first order the stability boundary is given by the 
condition 

Ac„ tt '(vo)dr 0 

This leads to the interpretation of the indicial curves shown 
in figure 3; namely', an airfoil pivoted at its leading edge can 




Oscillatory instability possible if 

c -;(oo) {jo)dT. 

Figure 3.— Interpretation of stngle-degree-of-freedoni stability boundary in terms of 
Indicial response cur res. 


- A pulse function occurs only for the case 3 /q= 0. Its treatment to expressions such as equations (5) fs disdussed to a subsequent section. 

4 That Is, the moment of Inertia fs large. Experimental verification of the existence of this type of instability can be found In reference 10. 



398 


REPORT 1077 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


have an oscillatory instability if the steady-state value of 
c m ' is less than the area between the steady-state value of 
c m J and its indicial value when expressed in chord lengths 
traveled. 

PART II— METHODS FOR SOLVING UNSTEADY LIFT 
PROBLEMS 

BASIC EQUATIONS AND SOLUTIONS 

In terms of the perturbation velocity potential <p, the wave 
equation can be written 

'Ptx J r<Pm J r , Pt*=-^i’Pt't' ( 8 ) 

Cto 

where a 0 is the speed of pressure propagation, t' is the time, 
and x,y,z are spatial coordinates. This equation applies to 
a flow field which is stationary at large distances from the 
disturbance region; furthermore, the coordinate system is 
stationary relative to the fluid infinitely distant so that, if 
a moving wing is being analyzed, the wing moves with 
respect to the x,y,z axes. The last equation can be put in 
a more convenient form by introducing the notation 


where M 0 =— <1, will satisfy both the above conditions. 

For example, suppose that the wing is moving along the x 
axis with velocity F 0 . Application of equations (10) to 
equation (9) makes the origin of the new system of axes also 
travel along the x axis with a velocity F 0 . This is seen to 
be so since £ is always zero when x=Vol'=M a t. Hence, 
the £ axis is “fixed” on the wing. As to the second condition, 
a straightforward exercise in partial differentiation yields 

%+^,+ % = ^rr (11) 

so that, in going from x,y,z to £,»j,f space, the wave equation 
remains invariant; consequently, both the requirements 
mentioned have been fulfilled. 

It is instructive to consider briefly the consequences of 
applying the Loren tz transformation. Although the wave 
equation remains invariant, such physical quantities as 
length and pressure do not. For example, a wing with a 
chord c in t he x,y,z space has, according to equations (10), 
a chord c/V 1 — A/ 0 2 in the £,ij,f space. Furthermore, the 
loading coefficient which, on the basis of linearized theory, 
is given in the x,y,z,t space by 


t=a4 r 


so that the dimension of t is length just as are the dimensions 
of the geometric variables x,y, and z. Equation (8), 
together with this transformation, yields the canonical form 
of the wave equation 

<Pxx+V’w+?’m = $ , «< (0) 

and it is this form which will be considered. 

The fust task is to study the relation between the motion 
of the wing and the coordinate system. As has already been 
mentioned, equation (9) is valid for a flow field produced by 
a wing moving relative to a coordinate system fixed with 
respect to the fluid at infinity. It is pertinent to consider 
the possibility of finding a transformation which will fix 
the origin of the coordinate system on the wing and, at the 
same time, retain the wave equation as the governing equa- 
tion of the flow. Certainly the first of these requirements is 
simple to fulfill if the second is neglected. 

The following transformation (known in relativity theory 
as a Lorentz transformation) 


i 




x — Mpt 

Vi=a? 


n=y 

r=z 

t—MoX 

r_ Vl~Afo s 


( 10 ) 


A p 4 dip 

qo^VoMo dt 

becomes for the £,ij,f,r space 

A p_^ 4 / 1 dip d<p\ 

2o = FoV1-A/o^ M o dr dt) 


( 12 ) 


(13) 


If the wing motion is steady and there are no transient effects, 
equations (11) and (13) are independent of time and, to- 
gether with the resulting length transformations, become 


%+^+%=0 

A p 4<P{ 

<Zo Fo a /f=17? 


Vi-M / V v ’ f Z J 


(14) 


These arc immediately recognized to bo Laplace's equation 
for incompressible flow and the familiar PrandtI-GIauert 
compressibility corrections. 5 

The preceding discussion has an important qualification, 
however, in the fact that the velocity of the moving co- 
ordinate system cannot exceed the speed of sound. A glance 
at equations (10) serves to verify this statement since those 
equations show M 0 must be less than 1 in order that £ and t 
be real for real x and t. In fact, it has been shown that there 
is no transformation which will fix the moving system of axes 
in a wing traveling at a uniform supersonic speed away from 
the original fixed axes and still keep the wave equation in- 
variant. Therefore, for analyzing a wing in supersonic 


‘ The equation r/V I— -W would first read i“(*— AfoM/Vl — Mtf where ft is a oonstant representing the ttme required for the motion to reach its steady state. However, the icoordl- 

Date can always be translated to any fixed position without affecting any of the equations for potential, loading, etc. Such a translation Is assumed to have boon made in equations (14), 
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flight, it is necessary to abandon one of the two proposed 
requirements: Either the coordinate system cannot be fixed 
in the wing, or the field equation must be modified. The 
latter of these two alternatives has been studied by several 
authors (see references 2, 3, 4, 5) but it is the former which 
will be considered in the present analysis. Further, since 
the axes cannot be made to travel as fast as the 'wing, they 
will not be made to move at all and equation (9) will be 
adopted throughout as the basic equation. 

Having decided upon the form of the partial differential 
equation, we must next establish the boundary conditions. 
For any given time these conditions are s imil ar to those 
studied in steady-state thin-airfoil problems; namely, either 
that the given slope of the wing surface is proportional to the 
vertical induced velocity <p, over the region occupied by the 
wing in the 2=0 plane, 8 or that the prescribed surface 
pressure is proportional to the timewise gradient <p t in 
velocity potential over the same region. The addition of 
time simply means that this region moves about in the 
2=0 plane in conformity with the direction and velocity of 
the wing. 

The solution to equation (9), subject to the boundary 
conditions just mentioned, can be expressed by a formula 
which may be regarded as requiring either the evaluation of 
a double integral or the solution of a double-integral equation, 
depending upon whether a boundary-value problem of first 
or second kind is considered. This solution is known as 
Ivirchhoff’s formula. (See reference 1.) It may be written 
in a form convenient for aerodynamic applications as follows: 


*= "w/J (t, 4 £r, 1 ‘- rJ - 

d yUi,yi.0,f- 




dS (15) 


where the A indicates the jump ( value on the upper surface 
minus value on the lower surface when applied to <p or j 


of the function in passing through the Zi=0 plane, 

r=y(ar— a-,)*+(j/ — yi)*-f-(z— 20 *, r„ =■/(*— »0*+(V— yi)*+s*» 
and where the area of integration S tt will be discussed in 
more detail later. 

The terms in the integrand of equation (15) can be 
shortened by introducing the following notation: 


(16a) 

d (dr d <p(x u yi,0,f— r)\ 

dz { r \dz,dr r ) 


dr n d TyT! 


(16b) 


In this notation, the subscript r in equation (16a) means 
that r is to be held constant in the differentiation, and the 
prefix A obviates the necessity of indicating that the func- 


tions considered are to be evaluated for z t =0, since it indi- 
cates that the difference of the values of the function across 
the 2 t =0 plane is to be taken. The right-hand side of 
equation (16a) can be recognized as a term representing a 
source located in the Zi=0 plane, and the right-hand side of 
equation (16b) is seen to represent a doublet located in and 
with axis normal to the 2 r =0 plane. The brackets [] about 
the functions in equations (16a) and (16b) have a special 
meaning which is defined in the following way: If / is a 
function the value of which at a fixed point P depends upon 
the coordinates x lt y u z u t of a moving point Q, so that 


then^ 


f=f(x l ,y l ,z l ,t) 


[f\=AxuVi,Zi,t—r) 


(17) 


where r is the distance from P to Q. As an example, con- 
sider the potential <p at a point P due to a moving source, 
the location of which at any time is <2- Then <p satisfies the 
condition just mentioned that it depends on the coordinates 
r,, Vu t of Q. The brackets [] indicate that the potential 
[<£>j depends not upon the source strength now at “time” f, 
but rather upon the source strength that existed “time” r 
ago. 7 For convenience, [<p] is referred to as the retarded 
value of ip. 

The expression for a doublet (equation (16b)) is usually 
expanded as follows: 


dr 

dz 




f (18) 


Finally, equation (15) becomes 




l^dr o 
r a dz 



dS 

(19) 


The application of equation (19) awaits only a discussion of 
the area S* over which the integration is to be made. This 
discussion is important enough, however, to merit considera- 
tion in some detail and will be given in the following section. 

THE ACOUSTIC PLAN FOBM 

Suppose that a line of sources is placed along the y t axis 
and that the strength of these sources is zero for f<0. At 
t = 0 they are “turned on” and, at the same time, start 
moving along the negative Xj, axis with the velocity T r 0 . 
After time V has passed, the source line has traveled a 
distance A / 0 t as shown in figure 4 (which is drawn for the 
case M o y> 1). Suppose next that there are two sensing 
elements, or detectors, placed at the point P(x,y ) located 
somewhere ahead of the y \ axis; one of these detectors is 
responsive to light and the other to sound. Now, the light 


• The i— 0 plane b assumed to be the “plane of the wtng ,r ; that Is, If the angle 0 / attack were lero and the wing had no thickness It would lie entirely In the i—0 plane* 

1 Quotes are used around the word time since the dimension of t Is actually length, not time. It Is convenient, however, to refer to f as “time," and, since the actual value of time la simply 
t divided by the constant at, thti should cause no confusion. 
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Numbers refer to t-r 



detector at any given time will show the sources lying in a 
straight line just as they would appear visually at each 
particular instant. The situation is entirely different, 
however, for the sound detector. First it is necessary to 
understand the nature of a spherical sound wave. Such a 
wave travels outward from its origin at a velocity do, so that 
in the time t' it has traveled a distance t. Before the wave 
reaches a point, the point is completely unaware (unaware 
is used in the sense that an instrument will record no change 
in any of the physical properties of the air at the point in 
question) of its existence and, further, after the wave has 
passed, the point remains subsequently unaware of its 
existence. Hence the only points disturbed by the wave 
are those momentarily on the spherical surface itself. (In 
this connection see reference 1, pp. 1-3.) 

The sound detector, therefore, can only “hear” sources 
which are so located that their spherical sound waves are 
just, at the given instant of time, reaching the detector. 
The locus of all the points which, at a time i ago, emitted 
sound waves that are just now reaching the point P{x,y) is 
itself a sphere and for convenience this sphere will be referred 
to as an “inverse sound wave.” 8 The traces of these inverse 
sound waves in the 2=0 plane are drawn in figure 4 as 
concentric circles about the point P(x,y). The intersection 
of an inverse sound wave of radius t — r with the line repre- 
senting the position of the sources at a time t'—r' ago gives 
the position of the sources which are just now signaling their 
presence to the sound detector at P. For example, when the 


source line started, it was lying along the y x axis. With 
reference to the present time t this was ayt r removed. Henco 
the intersection of the circle about P of radius t with the y x 
axis fixes the two points A and A', the spherical sound waves 
of which are now reaching P. A continuation of tliis process 
yields, for the locus of all points from which waves emanated 
that are just now toucldng the point P, the part of the 
ellipse shown in figure 4. 

In the sense that the light detector, because of the very 
large velocity of light, is “seeing” a straight line of sources, 
the sound detector, because of the relatively slow velocity 
of sound, is “hearing” an elliptic line of sources. Extending 
this concept to include a sheet of sources distributed over 
the surface of the wing, one can refer to the outline of that 
part of the wing which generates disturbances which can be 
measured by the light detector as the plan form (i. e., the 
visual plan form), and to the outlino of that part of the wing 
which affects the sound detector as the acoustic plan form. 
In a mathematical sense, the acoustic plan form is the area 
S a over which the integration of equation (19) is to bo made. 

The equation for the acoustic plan form can be formulated 
by means of the two equations 

(x-xry+iy-ytf+^it-ry ( 20 ) 

JiVu x u r) = 0 (21) 

where x, y, z are the coordinates of the point at which the 
induced effects are to be measured; x u y u the variable points 
of the sources; t, “time” now; and t—r, “time” ago. Equa- 
tion (20) is that of the inverse sound waves and equation 
(21) represents the position of the visual plan form at a 
“time” t. It is necessary to include the region behind the 
wing covered by the vortex wake as part of the visual plan 
form. In case the vorticity in the wake vanishes, as in the 
thickness problem, the wake may still be considered as part 
of the visual plan form, but the strength of the source-doublet 
distribution over that part of the acoustic plan form corre- 
sponding to the wake will vanish. If Q— 2*<0, the acoustic 
plan form does not exist since, for such a case, no source 
on the x u y x plane has had time to transmit its effect to the 
point x, y, z. On the other hand, if the circle (in the *i, y x 
plane) given by the equation 

(*— zi) s + (y— yi)*=^— 2 s 

lies entirely within the area occupied by the visual plan 
form of the wing at the beginning of the motion, the acoustic 
plan form is just this circle itself. For any other situation 
the acoustic plan form is formed in part, or in whole, by the 
curve found by eliminating r from equations (20) and (21) 
and in part, or not at all, by an arc of the circle {x— £i)*+ 
(y—y 1 ) t =t i —z i . One of the principal advantages of the 
acoustic-plan-form concept arises in the study of problems 
involving source or doublet distributions having constant 
strength. In such cases, the retarded values of the potential 


t The inverse sound wave baa for Its analogue In steady Ilf tlng-surfac8 theory the Mach foreoone. In that theory a disturbance outside the Mach forecone cannot affect the values of 
the induced velocities at the point where they are being measured. Similarly, in the present study, a source located outside the inverse sound wave of radius f cannot affect the values of 
any measurement made at the point P. • — 
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and its gradient appearing in equation (19) are constant 
and can be taken outside the integral signs. The problem 
is thereby reduced to the integration of a simple geometric 
variable over S a . 

A few examples will serve to fix the idea of the acoustic 
plan form. Consider first a two-dimensional, unswept wing 
moving at a constant supersonic speed in the negative r t 
direction. At time zero the leading edge of the wing was 
along the y t axis and now, at time equal to if, the wing has 
moved so that the leading edge coincides with the line 
Choose three points that are now lying on the 
wing. Let one point have its z coordinate in the range 
c z 2 >x >Vf*- z 1 (where c is the chord of the wing), 

the second in the interval y'f* — z*)>r> — y f®— z s , and the 
i bird in the range —\P—z*>x> — M 0 t. Designating these 
points by Pi,P t , and P t (see fig. 5), it can be shown that their 
acoustic plan forms are, respectively, a complete circle, 
a part circle and part ellipse, and a complete ellipse. The 
points P are at the centers of the circles and at focal points 
of the ellipses. Since, moreover, the circular plan form about 
Pi receives no signals from sources on the leading or trailing 
edge, conditions at Pi are consequently completely inde- 
pendent of the actual (visual) plan form of the wing. The 
elliptical plan form about P,, on the other hand, depends 
entirely on the shape of the leading edge; and finally the 
mixed plan form about P, is in certain regions (the circular 
portion) independent of the leading edge, and in other 
regions (the elliptic portion) entirely dependent upon it. 
Since the wing is traveling at supersonic speeds, the trailing 
edge and vortex wake can have no effect on the measure- 
ments taken on the wing and, in the same way, a point ahead 
of the wing leading edge, P 4 in figure 5, is undisturbed. 

Xext consider a wing moving at a constant subsonic speed 
in the negative X\ direction. As before, the leading edge 
was on the yi axis at t '= 0 and has traveled a distance — M a t. 
Choose now three points Pi, P 2 , and P a on the wing and 
unaffected by the wing tips. The acoustic plan forms for 
these points are combinations of circles and hyperbolas as 
contrasted with the circle-ellipse combination in the super- 
sonic case. Just as in the supersonic case, however, there 


. P . i 



Fificsi 5.— Acoustic plan forms lor points on wing traveling at supersonic speed. 


is a certain region represented by Pi in which the acoustic 
plan form is a complete circle and is independent of the 
visual shape of the wing (see fig. 6). Point P 2 is surrounded 
by a plan form which is part hyperbolic and part circular, 
the point itself being the center of the circle and the focus 
of the hyperbola. Point P a is a limiting value of P 2 ; 
it lies on the leading edge of the wing and the hyperbolic 
sides of its plan form have degenerated into straight lines. 
Finally, P 4 lies ahead of the wing; its plan form is still a 
combination of a hyperbola and a circle, but P 4 is now the 
focal point lying ahead of the hyperbolic branch used. 

Figure 6 was constructed so that the portion of the visual 
plan form behind the trailing edge had no effect on the 
potential at the various points Pi, etc. If these points had 
been chosen at positions where the wake could signal its 
effect, one of two acoustic configurations would result. 
First, if the wing is symmetric about the 2=0 plane, no lift 
is developed and the vorticity in the wake is zero so that the 
visual plan form need not include the wake, but effectively 
ends at the trading edge. In this case, the leading edge of 
the acoustic plan form is then determined as before, while 
its modified trailing edge may be made up, in part, of circular 
arcs formed by the primary wave and, in part, by an arc of 
the hyperbola formed by the (acoustic) intersection of the 
straight visual trading edge with the primary wave (such an 
arc being identical with the leading edge of the acoustic plan 
form but displaced backwards). On the other hand, if the 
wing has no thickness but is inclined to the free stream, it 
develops lift and the vorticity in the wake does not vanish; 
the acoustic plan form has a trading edge made up entirely 
of an arc of the primary inverse sound wave. The space 
between this arc and the acoustic trace of the visual trailing 
edge is covered by a sheet of doublets, the strength of which 
is determined by the vorticity distribution of the vortex 
wake. 

It is interesting to notice the conversion of terminology 
which arises in the analysis of unsteady lift problems. In 
the study of steady-state wings, it is customary (because of 
the nature of the governing partial differential equation) to 
speak of the subsonic problems as elliptic and the supersonic 



-.Wing at time zero 

Wing now 

FiacK 6.— Acoustic plan forms far points on wing traveling at subsonic speed. 
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problems as hyperbolic. Yet the acoustic plan forms just 
presented involved ellipses for the supersonic wing and 
hyperbolas for the subsonic case. 

To complete the remark, it can be observed that when the 
velocity of the wing is sonic the steady-state partial differ- 
ential equation becomes parabolic and, in this case, the 
acoustic plan form of a straight-edged wing also involves 
parabolas. When the leading edge is linear and normal to 
the stream direction, the eccentricity of the conic sections 
hounding the acoustic plan form for a point on the wing is 
equal to 1/M 0 . From this relation it is apparent that for 
M 0 less than, equal to, and greater than 1 the sections are, 
respectively, hyperbolas, parabolas, and ellipses. As might 
be presumed, the value of the eccentricity satisfies simple 
sweep theory so that, for an infinitely long straight leading 
edge, the eccentricity of the acoustic plan form is 1 /M 0 cos A 
where A is the angle of sweepback. The principal axes of 
the conic sections are always normal and parallel to straight 
leading edges. 


HOMOGENEOUS BOUNDARY- VALUE PROBLEMS 


Kirchlioff’8 solution to the wave equation can be applied 
to arbitrary wing plan forms undergoing arbitrary maneuvers. 
The boundary values for such general problems, however, 
usually lead to the development of double integral equations 
which are difficult to solve. As is usual in such cases, there 
are many special types of plan forms and maneuvers which 
lead to boundary-value problems that are simpler to analyze. 
An important class of these simplified problems is that 
arising from homogeneous boundary conditions. 


Let <p(x,y,z,t ) be a solution to equation (9). In certain 

special cases this can be written <p= (t) n <pa in which 

case <p is called a homogeneous function of degree n. The 
number of variables affecting <po is only three as compared 
to the four which are necessary to determine <p. If, there- 
fore, a partial differential equation can be set up for <po, it 
will contain one less mathematical “dimension” than the 


equation for <p. Following this observation it is necessary 
to proceed in two directions ; one to find the partial differential 


equation for <pq, and the other to find the physical problem 
and consequent boundary values leading to a homogeneous 
flow field. The latter path will be first explored. 

First, consider an example of a homogeneous boundary- 
value problem. Suppose that a rectangular flat plate starts 
suddenly from rest and moves forward at an angle of attack 
at a supersonic Mach number M a . At “time” the initial 


spherical wave generated by the forward right-hand comer 
has traveled outward to a radius f, and, at “time” 2t u to a 
radius 2t\. Figure 7 indicates the traces of these spheres in 
the 2=0 plane together with the original and present position 
of the wing leading edge. Let the points Pi and P a be 
located on the same rays through the origin of the circles 
and the wing comers. The problem is to find the pressures 
at P, and P a . 


It is apparent that, if every dimension in the figure involv- 
ing P 2 is divided by 2 £, and every dimension in the figure 
involving Pj is divided by t it the two figures will be similar 
in every respect and point Pi will coincide with point P a . 



Wing edge at time zero 

Wing edge at f, and 2t t 


Fioues 7 . — Geometric relationship for homogeneous flow. 

Since the vertical velocity iu 0 is constant over the plan form, 
a simple change in scale has made the boundary conditions 
for both problems identical. But this means that the solu- 
tions at Pj and P 2 are identical since the wave equation is 
invariant to change in scale. Hence, in regions of a rec- 
tangular wing unaffected by the waves from the trailing 
edge, the pressure can be written 

!<—«=!( Hu) <**> 

and the pressure is a homogeneous function of degree zero. 
A generalization of this example is contained in the following 
statement: 

(1) The pressure in any region affected by only two 
intersecting edges of a straight-sided flat plate traveling 
at a uniform subsonic or supersonic speed is homogene- 
ous and of degree zero (i. e., satisfies equation (22)). 
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Consider as another example the case of a flat rectangular 
wing traveling forward at a subsonic or supersonic speed and 
rolling about one edge, taken to be coincident with the x axis. 
The argument follows the same lines as before, and again a 
change in scale, proportional to the time, makes the geometr 3 r 
of the wing-wave combinations identical (in regions affected 
by only two intersecting edges) for different times. The 
boundary values over the wings will not be the same, how- 
ever, unless the slope of the w 0 distribution is adjusted in 
each case. But ir 0 can be adjusted by reducing it an amount 
proportional to the distance from the axis of rotation. The 
boundary-value problems are then similar for different values 
of time. Finally, therefore, the pressure can be written 


Ap , 

2o 



(23) 


which is a homogeneous function of degree one. A general- 
ization of tliis example is expressed as follows: 

(2) The pressure in any region affected by only two 
intersecting edges of a straight-sided flat plate traveling 
at a uniform subsonic or supersonic speed and rotating 
at a constant rate of pitch or roll is homogeneous and of 
degree one (i. e., satisfies equation (23)). 

It should be noted that both (1) and (2) are equally true 
for the steady-state case when all transient effects have 
disappeared. In supersonic wing theory they lead to conical 
and quasi-conieal flows, respectively while in the subsonic 
case they lead to flows about wings having infinite chordwise 
extent. In general, homogeneous flow occurs when the 
boundary conditions after a change in scale are proportional 
to their original values. 

Consider next the modification of the basic partial differ- 
ential equation (equation (9)) under the assumption that the 
flow is homogeneous. If the pressure is given by a function 
that is homogeneous and of degree zero, then, by equation 
(12), the velocity-potential function will be homogeneous 
and of degree one. If the notation 


^ £ 2/o> ^ Zo( 


a,yo, Zo)/ 


(24) 


is used, then equation (9) becomes 


( 1 -r„*)V 0 +(l -yo 5 ) $**>+(1 - zo 8 ) 

2j* g r 0 4 1 x 0 j 0 — 2y 0 Zo^ > ir 0 t 0 =0 (25) 

and a linear partial differential equation with three inde- 
pendent variables is therefore obtained. 

In the general theory of partial differentia] equations of 
second order, the character of an equation is determined from 
the geometric nature of a related quadric surface. The char- 
acter of equation (25) can be shown from such considerations 
to depend on the sign of the expression l—x^—y^—z^. It 
is immediately apparent, however, that within the unit 
sphere in the x^y^ space the sign of is 

everywhere positive and outside the sign is everywhere 


negative. It follows that outside the unit sphere equation 
(25) is hyperbolic and inside the unit sphere it is elliptic. 

The character of equation (25) is of particular interest 
since the difficulties inherent in the determination of the 
solutions can be estimated without actually obtaining the 
solutions. For example, consider the two configurations 
shown in figure 8. These wings started moving at f=0 with 
the foremost portion of their leading edges on the y Q axis and 
have by now traveled forward at a supersonic speed to attain 
the positions represented by the figure, the unit circle being 
in each case the trace of the primary wave from the vertex on 
the z=0 plane. Outside the unit sphere, the governing 
equation is hyperbolic and the behavior of the flow is similar 
to that in steady-state supersonic-wing problems. Inside 
the unit sphere, on the other hand, the character of equation 
(25) is elliptic. 

It is instructive to notice that this entire development has 
a direct analogue in the study of three-dimensional, steady- 
state, supersonic wings. In that case the original equation 
is the tliree-dimensional wave equation 

<Pxx—<Pn—<Pu = 0 (26) 


V 0 
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By considering the velocity potential to be homogeneous and 
of degree one, Busemann in reference 1 1 was able to introduce 
the transformations 




<?(x,y,z)=x$ (y a ,s 0 ) 


and transform equation (26) to the form 

(l-i/o*) ^Vo-Ki-Zo 2 ) $, <J , 0 -2y o2o ^ 0 , 0 =0 (27) 

which is the two-dimensional form of equation (25). Flows 
governed by equation (27) have become known as conical 
flows. A study of equation (27) shows it to be elliptic inside 
and hyperbolic outside the unit circle. In this case, however, 
the equation has only two independent variables so .that 
once the equation has been transformed (by means of the 
Tschapligin transformation) to the two-dimensional form of 
Laplace’s equation solutions are not difficult to find. 

The simplification of the four-dimensional theory brought 
out by the introduction of homogeneous flow was more ap- 
parent than real since the resulting partial differential equa- 
tion, although containing one less dimension, was unwieldy. 




(b) 




BOUNDARY-VALUE PROBLEMS INVOLVING NONINTERACTING SURFACES 


Another class of wing problems which is simplified both 
in theory and in practice by reasoning from physical knowl- 
edge of the flow behavior is that in which the wing has a 
supersonic edge (i. e., an edge which is traveling with a 
supersonic normal component of velocity). 

When the acoustic plan form is affected only by a super- 
sonic edge, it is apparent that the flow on the upper surface 
of the wing is independent of that on the lower surface. 
Hence the solution to such problems can always be written 
in terms only of sources as follows: 



where ^><p^>Z\=w u (xi,yi) = V 0 \ u (:Xi t yi) } X„ being the local slope 
of the surface in the direction of F 0 - Since the equation (28) 
is equally valid for symmetrical nonlifting surfaces and lifting 
plates, its value and simplicity are evident. 

If the wing plan form is further specialized by having not 
only supersonic leading edges, but also having a straight 
trailing edge perpendicular to the. direction of motion, addi- 
tional simplifications can be used.® Consider, for example, 
the two-dimensional wing (a) in figure 9. Let this wing have 
an angle of attack a(t) which varies with time in an arbitrary 
manner. There results from such an angle-of-attack varia- 
tion a certain lift which also varies with time. Hence, if 
L* represents the total lift on an airfoil of very high aspect 
ratio and c t * represents the section lift-curve slope, then 10 


L* 

2o<* 


=cc Ja *(c,()(span) 


rd-d{t) 



analysis is based on a linear partial differential equation, by 
superposition principles the total lift on wing (b) or (c) is 
just half of that on wing (a). In another sense, the lift 
coefficient for the whole whig based on the deflected area is 
the same for all three cases. A suitable superposition of 
wings (a), (b), and (c) will give whig (d), which then has 
the same lift coefficient based on deflected area. Finally, 
because of the supersonic stream, wing (e) can be obtained 
from (d), hence it also has the lift coefficient common to the 
other wings. It is, of course, necessary that the variation 
of a with the time be the same in each case. 

The preceding process can be extended one step farther 
to the development of the lift due to a single deflected 
element. By considering figure 10, it can be seen that 

1* c c, *{c,f)-(c—Ac)ci *(c—Ac,t) 

= = AS 

q 0 ct Ac 

where AS is the area of the deflected element and c is the 
distance Tram the centroid of AS to the trailing edge. By 
the usual limiting process the latter equation becomes 


Next it is clear by reason of symmetry that the total lift on 
wings (b) and (c) in the figure are equal. Then, since the 


U_ 

2o<2 


=^[cc (a *(c,f)]dS 


• The following method simply eitends, to include the effects of unsteady motion, a theorem given by Lagerstrom and Van Dyke. (See reference 12.) 
a The asteitsks on quantities indicate that two-dimensional values are taken, or that a hlgh-aspect-ratlo wing is considered and tip effects are neglected. 
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Fhii/EE 10.— Development of lifting element. 


Finally, if a wing is composed of a distribution of these 
elements, then, for a coordinate system centered at the 
apex of the wing leading edge, there results 

0 L = — ^ £ J a{ x, y) [(c a —x)c,*(c 0 — x, t)\dydx (29) 

where Co is the maximum chord (see fig. 11). In the develop- 
ment of equation (29) each element is assumed to have the 
same variation of motion with time. 


edges rising and sinking with a harmonic motion can be 
computed from a single integration of the results presented 
in reference 5. 

There is another simplified method for obtaining the total 
lift and moment on a wing with all supersonic edges and a 
straight tr ailing edge. In this presentation it will be assumed 
that the trailing edge is normal to the free stream. How- 
ever, since the wave equation is invariant to a rotation, it 
will be apparent that the solution can be generalized to in- 
clude a straight, supersonic trailing edge yawed with respect 
to the free stream. 

Consider equation (9) and integrate each term with respect 
to y between the limits minus and plus infinity. 11 There 
results the equation 




=0 



Notice that when the wing is a flat plate flying at a steady 
speed so that all transient effects have disappeared, c^* is 
independent of t and of the chord length, being, in fact, 
equal to 4//9. Then equation (29) becomes 

where a is the average angle of attack of the wing. This 
result has already been obtained in reference 12. When 
a(x , yi is independent of x (as for a flat wing sinking or 
rolling), equation (29) becomes 

b_ 

I" a(y) c c t *{c,t)dy (30) 


where c is the local chord which is, in general, a function of y. 
Equation (30) simply indicates that longitudinal strip theory 
is exact for calculating the total lift on such wings. 

Finally, notice that the calculation of the unsteady lift on 
three-dimensional wings with supersonic leading edges and 
straight trailing edges perpendicular to the free-stream 
direction has been reduced to an integration involving the 
relatively simple results for a two-dimensional wing under- 
going the same unsteady motion. For example, the lift 
on a flat unyawed triangular wing with supersonic leading 


If y=y l (x, z, t ) and y=y r (x,z,t ) are the equations of the Mach 
waves streaming back from the leading edges on the left and 
right sides of the wing, respectively (see fig. 12), then, since 
tp is continuous across these waves but <p x , <p„, <p t , and <p t 
are not, 


b ! 

bx 1 







where u r and u t are the values of u on the interior faces of 
the right and left Mach waves, respectively, and 



Values of the terms involving <p, and <p t are similar to those 
involving <p x so that finally, if 


then 


<$= 



br 2+ bs 2 bf 3 




(31) 


(32) 



Figcei 12.— Forward portion of Mach wave system for a su person le-edged, trfongnhr wing. 


it The basic idea for this solution was Riven by P. A. Lagerstrom. hi his lectures at the California Institute of Technology. 
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The terms enclosed within the parentheses in the last 
equation combine internally so that each is zero. For the 
case of interest here, this is not difficult to show. Consider, 
for example, the Mach wave streaming from the light edge. 
Then, since the equation of this wave is 

y r = — l+mx+mM^ 


and the value of the potential is the steady-state two- 
dimensional value given by the expression 


m 


<P,— — w u - 


(x+M a t)— y — 1 




the term 


becomes 


and this is identically zero. 1 * Finally, therefore, equation (9) 
has been reduced in terms of equation (31) to 

(33) 


The boundary condition in terms of $ for a triangular 
wing with supersonic edges is given by 



where w u is the vertical induced velocity in the plane of the 
triangular wing. The derivative with respect to z can be 
carried tlirough the integral sign because the extra terms in- 
volving the value of <p at y\ and y, vanish. In fact, since 
equation (34) applies to the 2=0 plane, the limits y ( and y r 
can be replaced with the expressions for the left and right 
leading edges of the triangular wing, respectively. The 
boundary condition expressed by equation (34), used in con- 
junction with equation (33), suggests a problem exactly like 
those posed by lifting surfaces in steady-state wing theory; 
in fact, the problem of a wing tip of specified camber in a free 
stream at Mach number ^2. Figure 13 shows a lifting sur- 
face in the x,i plane. The solution for the potential $ in the 
steady-state problem can be written 




1 


x 



b* 


bz 

t-0 

VG-ti y- 

-(x— Z,) J 


dtidx i 


where a is the area on the wing plan form in figure 13 that is 
included in the forecone (f— #i) 2 = (a: — Zi) 2 . 

Now from equation (12) we have 



4 

VoM 0 bt\ 


, — 0 



Fiqure 13.— Region of xt plane In which boundary conditions for* aro known. 


whore s is the local semispan of the triangular wing. This 
equation shows how the solution to the lifting surface prob- 
lem in * will aid in the solution of the unsteady problem, for 
the unsteady lift on the triangular wing is given by 



4 

SM 0 V a 


f c ° bd> 

Jo & 


dx 

i-O 


where S is the area of the triangular wing. It is therefore 

^(J> , 

seen to be convenient to evaluate the quantity , given 
by 

a... 


* — 0 


bt 


Jl 5 

X 


4 f f dt,dx x (3 5) 


The pitching-moment coefficient for the triangular wing in 
unsteady motion can be evaluated similarly. Specific appli- 
cations of this method will be found in a subsequent section. 


TWO-DIMENSIONAL BOUNDAET- VALUE PROBLEMS 

The simplification brought about when the flow is inde- 
pendent of one dimension is obvious. In such cases, the 
three-dimensional wave equation (9) reduces immediately to 
the two-dimensional wave equation. Typical examples of 
this type of problem can be constructed by considering flat 
plates wdiich start suddenly at t—0 and travel thereafter at 
constant supersonic velocities. Two examples, one a corner 
of a rectangular wing and the other a triangular wing, are 
showm in figure 14. After time f=0, the edges of the wings 
send out cylindrical waves and the outer boundaries of these 
waves at time t are shown as dashed lines parallel to the edges 
in question. Since points in regions 1 and 2 are affected only 
by a single edge, the wave phenomena in these regions are 
cylindrical, and the physical quantities are in both cases 
independent of distance parallel to the edge which acts as 
their generator. Hence, the flow field in these regions is two- 
dimensional. (Region 3, incidentally, is independent of 
distance in both x and y directions and is, therefore, one- 
dimensional.) 


i* It is not necessary to perform a direct calculation In order to prove the above result for arbitrary plan forms. The terms in parentheses In equation (32) represent the directional deriva- 
tive of the velocity potential tafcen along the so-called "oonormnl” of the foremost disturbance surface. Since e Is constant on the surface, and since the cunormal lies along the surface, the 
terms in parentheses are zero. 
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Solutions to the two-dimensional unsteady problems are 
sometimes especially easy to find because of the analogy they 
have with three-dimensional, steady-state, lifting-surface 
problems. For example, consider an infinitely long unyawed 
wing which starts from rest and travels forward at a velocity 
V 0 which may or inaj- not be a function of time. The trace 
of this wing in the x, t plane is like that shown in figure 15. 
(In the figure shown, the wing velocity is varying and is 
always less than the speed of sound.) The boundary condi- 
tions are that <p, is specified over the shaded area and the 
loading Ap, is zero everywhere except within the shaded 
area. But if x is replaced by y and t by x, these boundary 
conditions are exactly the same as those for a plate of known 
camber and angle of attack, with a plan form as indicated 
by the shaded area, placed in a free stream directed along 
the positive x axis at a Mach number equal to \'2. The 
solution for the one problem may be used, therefore, as a 
solution to the other with only a change in notation. 


- Traces of characteristic cones 



BOUNDARY CONDITIONS FOR VERY SLENDER WINGS 

When the wing plan form is slender in the sense that its 
length in the streamwise direction is large compared to its 
span, an estimation of the loading on it can be obtained by 
neglecting in the partial differential equation the gradient 
of the induced velocity component in the stream direction. 
Thus, if the wing is moving in the negative x direction, 
equation (9) reduces to 

(36) 

which is again the wave equation but in two space dimen- 
sions. Since equation (36) is independent of x, study can be 
made independently in each plane x= constant. This is an 
extension of steady-state slender wing theory, see, e. g., 
reference 13. Figure 16 (a) shows a typical section in the yt 
plane. If the wing is a flat plate at a constant angle of attack, 




(a) Om fixed relative to 5tUI air at Infini ty. 

(b) Axes fixed on wing. 

Figcbe IS.— Slender win* In unsteady flow. 


Fiui’&i IS. — Decelerating wing in it plane. 
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the value of <p, over the part of the yt plane occupied by the 
wing is a constant and the jump in <p across the vortex wake 
must be consistent with its value at the wing trailing edge. 
The analogy with three-dimensional, steady-state, super- 
sonic, lifting surface theory is apparent. 

There is another way by means of which the effects of 
unsteady motion on slender wing pressures can be estimated. 
If instead of using the stationary xyzt coordinate system, the 
reference axes are fixed on the wing by the simple set of 
transformations 

Xi=x+M 0 t 

ti=t 

2/i=l/ 

Z\ = Z - - 

equation (9) becomes 

( 1 —Mo*) <«*,*, — 2M 0 ip Zx x t + <py lVl 

Again, if the induced velocity components in the stream 
direction are neglected, the simplified equation 

< Py l v l J r’Pt l t l =<Pt 1 t l 

results. The latter equation is identical in form to equation 
(36). However, now the axes are fixed on the wing and a 
typical section in the j/ifi plane is similar to that shown in 
figure 16 (b). In this case, a flat plate wing is represented by 
a constant value of over the entire shaded area in the figure 
and elsewhere <p must be continuous. 13 

TWO-DIMENSIONAL UNSTEADY INCOMPRESSIBLE FLOW 

The analogy between two-dimensional unsteady and three- 
dimensional steady flow includes the case of a two-dimen- 
sional, unsteady, incompressible flow field the analog of 
which is a three-dimensional, steady flow field having a 
free-stream Mach number equal to 1. Tliis can be demon- 
strated by inspecting equation (8). Since the flow is two- 
dimensional and since for an incompressible medium the 
speed of sound a 0 is infinite, the basic equation governing the 
flow can be written 

<Pxx-\-<Ptt = 0 

It must be remembered, however, that time still appears in 
the boundary conditions and in the equation for the loading 
coefficient which, according to equation (12), can be written 

A p 4 b<p 

?o Vo* bt' 

Hence, the basic partial differential equation and the expres- 
sion for the loading coefficient are the same, except for a 
change in notation, as those governing three-dimensional, 
steady-state problems when M 0 =l. 

If a two-dimensional wing in an incompressible fluid 
starts from rest and travels forward at a speed V a> the trace 
of the wing is as shown in figure 17. The essential difference 



between this problem and the more general case of two- 
dimensional compressible flow lies in the fact that in this 
case the traces of the characteristic cones are normal to the 
t' axis. The boundary conditions are therefore satisfied 
along lateral strips and, in lifting-surface terminology, the 
analysis corresponds to slender-wing theory. These latter 
methods are well established and in reference 13, for example, 
the manner in which the Kutta condition is imposed is dis- 
cussed in some detail. The trailing vortex sheet for the 
lifting wing has the same distribution of vorticity that exists 
behind the unsteady two-dimensional airfoil and the rolling 
up of the vortex sheet can be studied from either standpoint. 

PART III— SOME APPLICATIONS OF THE METHOD 

STARTING LIFT OF A WING 

One. of the simplest and yet most general respite which 
can be derived on the basis of the present theory is the initial 
value of pressure on a wing surface starting suddenly from 
rest with a velocity T r 0 . The discussion will be made for a 
wing without thickness although the method will be seen 
to apply to the thickness case as w r ell. 

Consider a surface with a plan form as indicated in figure 
18. The acoustic plan form of a point P(x,y) on ihe surface 
is a small circle of radius t. Since no point on the wing out- 
side this circle can influence the pressure at P, the upper 
surface is independent of the lower surface, except for a 
band of width t around the edge of the wing. It is, there- 
fore, evident that the boundary-value problem to be solved 
has been treated in the section Boundary-Value Problems 
Involving Noninteracting Surfaces. The solution follows 
directly from equation (28) and can be written 



u The analogous problem In steady-state wing theory is that of a low-aspect-ratio, rectangular, flat plat* in a free stream having a Maoh number eaual to 
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Wing plan form at t=o 
Wing plan farm now 



Using a polar coordinate system defined' by 

x—x x —r cos 0 
y—yi=r sin 8 
dxi dy t =rdrd8 

there results 14 

<P** — Wu !f’ y) f d8 f dr=—tw a (x,y) 

JO ,/Q 

so that 

(±P\ |_ i /V\ _ 4 «? u (j, y) 

V ?oA-o t T 7 0 MoU/,.r VoMo 

If a=—w x (x,y)IV a is the local slope of the wing, the expres- 
sion for load coefficient becomes 


/Ap\ _ 4a(x,y) 

\9oA-o A/q 

The starting value of lift coefficient can, therefore, be 
written 



where a is tbe average angle of attack of the surface defined by 



(S’ being the area of the wing plan form. 


INDICIAL FUNCTIONS FOE A TWO-DIMENSIONAL FLAT PLATE, Mi-9 

The basic partial differential equation governing an 
unsteady flow field in the two-dimensional incompressible 
case is Laplace’s equation in two dimensions 


0 (37) 

where x is distance along the chord and z is height above the 
plane of the wing. It was also pointed out in Part I that 
the boundary-value problems arising in the study of the 
unsteady, two-dimensional plate flying in an incompressible 
medium were directly analogous to those which are studied 
in three-dimensional, steady-state, lifting-surface theory 
under the classification “slender- wing theory.” This analogy 
is useful since well-established concepts in one field can be 
immediately carried over into the other. It should be 
mentioned, however, that the subsequent treatment of the 
incompressible case is not intended to be an improvement 
on Wagner’s original derivation (see reference 14) but 
rather it is a rederivation along lines that will be used later 
in the analysis of the compressible case. 

The initial pulse —The first analogy with slender-wing 
theory which will be used concerns the initial pulses that 
occur in the values of lift and pitching moment. It is a 
well-known result (reference 15) that the total lift, as given 
by slender-wing theory, on the wing shown in figure 19 (a) 
is a function only of the maximum span and the value of 
Wu along the section of maximum span (section AA). It is, 
therefore, independent of the wing twist and leading-edge 
shape ahead of section AA. This concept has been extended 
in slender-wing theory to the extreme case shown in figure 
19 (b) of a rectangular wing. The lift on such a wing is 
concentrated entirely along the leading edge and is a function 
only of the span of that edge and the value of there. By 
the analogy existing between the two theories, therefore, it 
is evident that the solution to the indicial problems in two- 
dimensional, incompressible, unsteady flow (fig. 1& (c)) will 
contain a pulse at f'=0. 

The evaluation of this pulse will be treated briefly. A 
solution to equation (37) for the vertical induced velocity 
in the 2=0 plane can be written in terms of the jump in u 
across the 2=0 plane (see reference 13) ; thus, for the shaded 
area in figure 20 this is 

1=7; dx ' <38) 


The general inversion of equation (38) can be written 

A 


A «(x) 


where 


«V(x-|-a)(6 — x) 


■f 


nr^{x-\-a) (b — x) 


= r “^V(*i+a)(*-*0 dx, (39) 

— X) * — x l 




A udx 


» The mean value theorem gives 


(£,q) 

2x 



vrhere £ and * lie somewhere in 


Hence, m f approaches 0 t £ and * approach c and jr r respectively. 


27*2483—54 27 
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U 

(b) 


Pulse m lift 
at t'0 — ' 

\ 



Fiooks 19.— Lift pulses In slender- wing theory. 


In the present case A is zero, since A <p is zero at x=—a and 
i=b, and an integration of both, sides of equation (39) with 
respect to x between the limits b and x gives 


A<p(t', x)=—— f w(t',x 0 

irj-a ; 

rn y’(6-a' 1 )(g+z)+N/(a+y 1 )(&-j:) dx 
s f (x—xi )(a+6) 1 

Adoption of the notation 


v{t', z 1 )+constant= 


J* w(f xi)dxi 


and integration by parts leads to the equation (since 

f vanishes if &>z> — a) 

J-afr-XiHib—xJia+X!) 


A <p= 


2rJ(b—x)<ja+x) 


£ 


v(t', xt)dxi 


(x— Zi)V(6 — xOfa+zO 
(40) 


The loading can now be determined by using equation (12) 
(and differentiating with respect to t' rather than t) 

Ap_ 4 £>e> 

to~v a *w 


If the shaded area in figure 20 is allowed to vanish, all the 
loading accumulates along the x axis in the region 0<z<c. 
Therefore, the integral of the loading with respect to V over 
the shaded portion must be considered. The final result for 
the pulse loading (Ap/q a ) s at f'=0 can be expressed in terms 
of the 5 function (see list of symbols) as 




1 y( 0, *i)dxi 

(Z — Zi)V(c— Xi)Xi 


(41) 


The boundary conditions for the sinking and pitching wing 
given by equations (2) and (3), when inserted into equation 
(41), yield 


48(0 

Fa 


V(c— x)x 


m = 

\g.oa/t 

1% 


\qodJi 


F« : 


(42) 


After integration, the pulse values for lift and pitching 
moment may also be obtained. Hence, 


( C '“) i= 2Ti Kn 

(•O.— 4 Tl^) 

( C " e ) #= ^?4 


(43) 


where the primes indicate that the wings are pitching about 
and the moments are measured about the leading edge. 
These expressions may be inserted in equations (4) and, since 
the integrand becomes zero everywhere except at the point 
the contributions of the initial pulse to the expressions 
for the lift and pitching-moment coefficient developed by an 
arbitrary variation of a and 6 with time are 


ire . . xc s 
Ct »~W 0 a+ IVo* 


9 


, ire . 9_ tc 1 •• 

Cn >- Wo a 04 TV 

where a and 9 are evaluated at t'. 



Fig u ns 20 .— Undegenerate form of Initial pulso 
In lift. 
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The variation for — The integral equation (38) is still 

perfectly valid when applied to the flow field for £'>0. 
It is convenient to rewrite the equation in this case, however, 
so that the effects of the vorticity on the wing and in the 
wake are separated. Thus, 


w u = 


J_ *'0*' x t ) 1_ r c Au*(Xj) 

2 xJ-r 0 t' a: — jti 1 27r,]c-v 0 £' x — x t 


dii 


(.44) 


where A u* (xi) is the value of Au in the starting vortex wake. 
It is independent of l' since its value at all points along the 
line ab in figure 21 is the same as at the point a. 



Since the value of A < p is zero at the leading edge and at the 
trailing edge is equal to the total circulation r, two altema- 
tive forms for the lift can be written 

• - d Ce-V 0 t' 

l=PoV<F+p a jjj j_ V ' tt A *P,x)dx (48) 

and 

Jr A rc—V^t' 

i=P*(c- WO |_i V . *A«(f ',x)dx (49) 

By substituting equation (45) into (49) and integrating, it 
can be shown that 


l=Tp a cVc, 2 Cf 


PqcVq r e Au* ( xi ) 

2 Jc-Vot' V(x x + I VO (Xi-e+ WO 


dx i 
(50) 


Since the value of A«*(xi) has been determined by Wagner, 
equation (50) can be evaluated to obtain a solution for the 
section lift. A plot of the section lift coefficient is shown in 
figure 22. Initially there is the pulse having an intensity 
defined by equations (43). After the pulse at £'=0, the 
value of the section lift coefficient starts at one-half its 
asymptotic value. It then increases, slowly approaching its 
asymptote of 2-xct. 


A reduction of equation (44) can be obtained for the case - 
of the sinking wing, where w x — — V Q a, by using the inversion 
given by equation (39). Thus, 

x Au(t r , x)y(x + VoO(c— Vtf— x)=A— jr(2x+2Fo£'— c)F 0 q:-|- 



, V(*I+ W0(*i-«+ VST] 


1 (x-xO J 


Since A is given by the relation 

A= J ^ (/ Au(f,xj)dxi=— Au*(ar,)ixi 

i 

it follows that 


w „ T>(2 j+2F/-c) , 1 

V(*+W0(c-W'-Jf) ^(x+V^c-W-x) 

P VCxi+TWi-c + TO dx, (45) 


According to the Kutta condition Au(t'^e) remains finite as x 
approaches c—V<t'; the integral equation for Au*(x,) thus 
becomes 


h £*• VYvTT? ix ' (46 > 


which was derived and studied originally by Wagner (refer- 
ence 14). 

The section lift and pitching moment can be derived in 
terms of Au*(x t ) in the following manner. By definition, 
the section lift l is 



By definition, the pitching moment can be written 

t=-| poTY r V dV+r) ^-dx 
2 J-v 0 t' q 9 


m= 


rc-Vot’ 

=— /»oJ_ tv , (V 0 t'+x)^-dx (51) 


where the moment is taken about the leading edge and is 
positive when the trailing edge is forced down. A develop- 
ment, similar to the one given for the lift, gives 


, 1 i-i A P / f* -1 * 1 8 Ap , .... 

4=5 Pol o 1 — dx = p a -5-r dx (At) 

2 J-ioi' q a J-Vtt' be 1 



( 52 ) 
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This result, that for £>0 the indicial center of pressure re- 
mains constant at the quarter chord throughout the motion, 
is classical. 

If the boundary condition for a pitching wing, .w u = — 
(ar+ V a t') 6, is substituted into equation (44) and the inversion 
given by equation (39) is again used, it can be shown in the 
same manner used in the derivation of equation (46) that for 
x=c— V 0 t' the relation 


3 

4 


r au * 

njc-Vot' 



gi + V 0 t' 

Xj — C+ Vat 


idxi 


By means of equations (52) and (57) the expressions of 
the three iudicial functions, c m J, and c n can all be 
written in terms of the indicial lift function for t>0. Hence, 



(58) 


applies. This integral equation applies to a wing pitching 
about its leading edge. If, instead, the wing is pitching 
about the three-quarter-chord position, an essential simpli- 
fication is achieved. In this latter case, downwash is given 
by the expression 

w u = — 6^x+Vot'~ (53) 


The variations of the four indicial functions will be shown 
later in Part IV. For values of r 0 (,= Vat'/c) larger than those 
shown in figure 22 the approximate equation suggested in 
reference 16 can be used, namely', 


and the resulting integral equation becomes 

0= P Au,* (?,) J -ll+Vg- d Xl (54) 

,Jc-v„r ^ V «i — c+Fot' 

where Ati**(a!i) represents the vorticity in the wake following 
such a motion. The solution to equation (54) is simply 

A«,*(i,)=0 (55 

From equation (55) it follows that the total indicial lift for 
£'>0 on a wing pitching about the three-quarter chord point 
is zero, and that the wing wake is free of vorticity. Further, 
it can be shown that the total indicial pitching moment (still 
measured about the leading edge) is 


This alternative result has, according to reference 16, an 
error of 2 percent or less for the entire range of time from 
0+ to infinity. 

INDICIAL FUNCTIONS FOB A TWO-DIMENSIONAL FLAT PLATE, A/t-O.S, 0.1! 

When the Mach number is no longer small, the analysis 
in the preceding section must be modified. As an example 
of this modification, we shall evaluate the indicial response 
on a sinking wing flying at a Mach number equal to 0.5 and 
the indicial response on a sinkmg or pitching wing flying at a 
Mach number equal to 0.8. 

Since the wing is two-dimensional, the partial differential 
equation governing the flow field (equation (9)) reduces to 

<Pa+<Pi,=<Pit (60) 


m 


pice 3 

TiT 


V 0 6 


(56) 


The transfer of equations (55) and (56) back to the case 
in which the wing is pitching about its leading edge can be 
readily' accomplished by means of the boundary condition 
shown in figure 23. Hence, if (c {j ) a/ 4 refers to the lift coeffi- 
cient on a wing pitching about the three-quarter-chord point 
and 4 refers to the pitching-moment coefficient meas- 

ured about the leading edge of a wing pitching about the 
three-quarter-chord point, then 



where it must be remembered that the axes are fixed with 
reference to the still air at infinity and the wing is moving 
in the 2=0 plane. The equation for the loading coefficient 
remains as in equation (12). The analogy which existed in 
the incompressible case between the theory for the unsteady, 
two-dimensional wing and slender-w'ing theory exists in this 
case between the theory' for the uusteady two-dimensional 
wing and the theory for a steady-state, three-dimensional 
wing traveling at a supersonic speed. Thus, in the three- 
dimensional, steady-state case the partial differential equa- 
tion governing the flow is 


<PnV'P,i — $ i <{>xx 


(61) 


and the equation for the loading coefficient is 


Ap_ 4 bip 
q 0 V 0 i>x 


(62) 


The boundary' conditions are in both cases that <p, is given 
over a portion of the 2=0 plane. It is evident by a compari- 
son of equations (60) and (61) and equations (12) and (62) 
that results from the three-dimensional, supersonic, steady'- 
state case (hereinafter referred to as the steady-state case) 
can be transferred to the two-dimensional, unsteady' case 
(hereinafter referred to as the unsteady' case) simply by 


Figure 23.— Change in boundary conditions corresponding to cl Lange in pitching axis. 
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replacing x, y, and jS in the former case by t, x, and 1, re- 
spectivelv, and by dividing the result for the loading coefficient 
by d/ 0 . 

The analog to the boundary condition for the problem of 
finding the indicial loading on a two-dimensional wing flying 
at a subsonic Mach number (fig. 24(a)) is the boundary 
condition for the problem of finding the loading on a constant- 
chord, swept-forward wing tip with a subsonic trailing edge 
such as that shown in figure 24(b). The Mach cones in the 



(b) 

i U nsteady case, 
fb) Steady-state case. 

Fir, r re 24,— Boundary conditions for a two-dimensional unsteady wins: mov- 
ing at subsonic speed and the analogous three-dimensional, steady-state 
wing. 


steady-state case, traces of which are shown as dotted lines 
in figure 24(b), become, in the unsteady-state analog, the 
locus of the sound waves which started at f=0 from the 
leading and trailing edges of the two-dimensional wing 
(fig. 24(a)). Finally, the analog in the steady-state fi eld 
of the unsteady wing would be a flat plate for the unsteady 
sinking wing and a plate with a linear variation of twist for. 
the unsteady pitching wing. 

Just as in the section on incompressible flow, the analysis 
will be divided into two parts. In cases for which J/o^O, 
however, the indicial functions contain no pulse at i=0. 
Hence, the first part of the study will be concerned with 
the behavior of the indicial functions in an interval for 
which t is small but finite and the second part, with their 
asymptotic behavior. 

The early stage. — The analog which exists between the 
steady-state and unsteady cases may be utilized to great 
advantage since the special methods and techniques developed 
for the solution of problems in the former case may be applied 
to the solution of the analogous problems in the latter 
field. In this manner an exact solution for the loading over 
the first five regions shown in figure 25 was obtained for a 
Mach number equal to 0.8, and for all the regions indicated 
for a Mach number equal to 0.5. Solutions for larger 
values of the time could also be calculated, but the labor 
involved in computing such cases becomes prohibitive and, 
as will be shown later, approximate methods can be developed 
which extend the solutions for the indicial lift and pitching- 
moment curves to their asymptotic values. 



Figi. ei 25.— Regions used In analysis of subsonic unsteady win*. 


The analysis used to calculate the loading over the regions 
shown in figure 25 is outlined in appendix A. Plots of the 
indicial loading on a sinking and a pitching plate flying at a 
Mach number equal to 0.8 are shown in figure 26. At 
time equal to zero the loading is constant for the sinking 
wing, and as time increases the loading-coefficient curve 
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(a) Sinking wing. . ' 

(b) Wing pitching about Its leading edge. 

Siouee 29.— Variation of two-dimensional Indldal load distribution with percent chord 
lor a Mach number equal to 0.8. 

it A similar result was noted In the stud; o Ithe load distribution on swept-back wings 


approaches the familiar two-dimensional, steady-state shape 
given by the equation 

Ap_4a I c — x 
Qo~ 

where x is distance from the leading edge which is at u. 

Figure 27 shows the load distributions for 4fo=0.5' and 0.8 
at the last value of time for which the exact loading curve 
was calculated. Notice that in each case the distribution is 
essentially the same as that obtained at time equal to infinity 
(i. e., the agreement is good with the curve produced by 
multiplying the right side of equation (63) by a constant 
factor). 15 .“The use of this fact simplifies the subsequent 
analysis concerning the asymptotic behavior of the indicial 
curve, i 1 

The indicial lift and pitching-moment functions were also 
calculated (see. appendix A) up to the time r 0 =2.333 for the 
wing flying at i/ 0 =0.5 and r 0 — 4 for the wing flying at 
A/ 0 =0.8. Their variation in this interval will.be shown hi a 
subsequent figure (fig. 59). It is evident from a glance at 



(a) Aft-OJS. 

CbJ Afo-OS. 

Fiouee 27.— Load distribution at end of early stage. 


subsonic leading edges (reference 17). 
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the figure that the calculations must be extended beyond 
this early stage since the asymptotic values are not even 
closely approached. 

Before studying the nature of these curves for large values 
of time, however, it is useful to examine them with reference 
to the discussion in the previous section on incompressible 
flow. For example, it was pointed out that the indieial 
center of pressure on the sinking wing remained at the 
quarter-chord point for time greater" than zero. Hence, 
let us consider the location of the center of pressure on the 
sinking wing when the Mach number lies between 0 and 1. 
By means of the indicia! curves for and c m J and by the 
relationship 

= (r/c) c .».Cj a 

the variation of (xlc) c . r . is easily evaluated ((x(c) e . p . is the 
distance between the leading edge and the center of pressure 
divided by the total wing chord). This variation is shown 
for the two Mach numbers in figure 28. It is apparent that 
the center of pressure is very close to the quarter-chord 
position for values of time greater than those for which the 
exact calculations were carried out. In other words the 
significant effect of compressibility' on the location of the 
center of pressure is limited to the interval ,0<t a <2 for 
d/„=0.5 and to the interval 0<f„<5 for A/ a =0.8. 

Likewise, it is apparent from the discussion of the incom- 
pressible case that the indieial functions for the pitching 
wing can also be expressed in a more convenient form by 
shifting the axis of rotation from the leading edge to the 
three-quarter-chord point. The values of and 

f°r 11 / 0 = 0.8 were calculated from the definitions 
given in equations (57) and are shown in figure 29. In- 
spection of this figure again shows that at a Mach number 
equal to 0.8 the compressibility effects are limited to the 
interval 0<f a <5. 



Knr-ai 28.— Center-of-pressure variation on sinking wing during early stage. 

The later stage. — It follows from the preceding discussion 
that when is large, the values of the indieial functions 
c„ a ', Ci', and c^' for compressible flow can be expressed in 
terms of c ta by equations similar to equations (58) which 
were derived for incompressible flow. Thus, after several 




Figure 20.— Variation of lift and moment on pitching wing during early stage; .V<— 0-S 


chord lengths have been traveled, one can write, on the 
basis of the asymptotic values shown in figures 28 and 29, 


Cm a '=*-cji 

c */ =3c *J 4 

c'=-(Scjie)-(ir!8p) 


(64) 


It remains, therefore, to determine the asymptotic behavior 

Of Cf . 

a s - 

Consider the steady-state solution for the lift on a two- 
dimensional, Sat lifting surface traveling at a subsonic 
Mach number. As was pointed out by 'W'ieghardt (reference 
18), if the lift on such a surface is represented by placing at 
the quarter-chord point a vortex which has the same circu- 
lation as that developed by the wing, the angle of attack 
measured at the three-quarter-chord point will be the same 
as that of the fiat plate. Extending this concept to include 
the unsteady effects, an investigation will be made of the 
variation with time of the vortex strength which will main- 
tain a constant angle of attack at the three-quarter-chord 
station following an impulsive start at i'=0. 

The analogous problem in steady-state theory becomes one 
of fin ding the strength of the vortex system, shown in figure 
30, which gives a constant value of w along the line CD. 
In the vicinity of the origin, of course, this representation 
gives a poor approximation to the original boundary-value 
problem. On the other hand, in this vicinity the exact values 
of lift and moment have already been determined. 

Each vortex composing this system lies along the line AB, 
extending from minus infinity toward the origin, and trails 
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back parallel to the x, axis to form the trailing vortex sheet. 
Note that, for convenience, the origin of the axis system has 
been located at the quarter-chord point. The solution to 
such a problem in steady-state, lifting-line theory would re- 
sult from the solution of the integral equation 16 (for a de- 
velopment using the notation adopted here, see reference 19) : 



AyiyJdXj 

[(x-xtf y l ) 1 \ itl 


where A<p is not a function of x L since the strength of a trail- 
ing vortex is, of course, constant. The area of integration 
r is the region within the forecone springing from the point 

x,V- 

If the above equation is transformed by means of the 
analogy to represent the solution of the unsteady problem 
(see fig. 30(a)), /S, x, and y are replaced by 1, f, and x, re- 
spectively; and Aip, the total jump in potential at a given 
section, is replaced by the circulation r. Hence, 




(a) Unsteady case. 

(b) Steady-state analog. 

Ftiii.tRs 30. — Vortex systems for two-dimensional unsteady and analogous 
three-dimensional, steady-state wings. 



rtr,) 


dt l 


where r, as indicated in figure 30 (a), is the area in the fore- 
cono from the point P which lies always along the line 
x=(lJ2)c—M 0 t. Integration with respect to U reduces the 
last equation to 



which, by means of the substitution x l jc=—x i j2 becomes 
c 

along the line 

_.l f x o (Xo+iUQ-galrCra) 

10 crpj 0 Ois-Xo+^VlXo-^CXo+Mi-rj )^ 2 C05) 

where \g=2Al 0 tg — /r 0 , n 0 =A[ 0 l(l-\-A{ 0 ), ^ l =2il/ 0 /(l — A/ q 2 ), 
/U}=l/(l+il/ 0 ), tg=t/c, and where, of course, w is a constant 
equal to — Tott. 

A solution for r(j 2 ) in the integral equation (65) may be 
obtained by expanding T in a series of the form 


rw 


irCaVg I ^Xj 

0 L V (ao+2j) 




x 2 


(ai + Xif 



Place equation (66) into (65) and expand in powers of 1/Xo. 
There results the expression 


in which 


w . Ci . c^n(l/X q) Cj 
IV r X 0 ^ Xo 2 ‘ X 0 * 


Cj = f>i 


£o_j 1_ 

2 ' /3* 


(07) 


Hence if a 0 and are chosen so that Ci is zero, an expression 
for T will be obtained which represents the solution to the 
integral equation (65) correct to the first order in l/t 0 (i. e., 
1/Xo) for large values of tg. Further, if equation (CG) is 
expanded in powers of j 2 , there results 


r W -«2l ! [, + L( 6l -|)+...] 


c c 

*• The symbols j* and j- arc used to Indicate that the finite part Is to bo taken. 

ni)d f _ „ a r • Art . 

“VP* a?J. 


nice Integral of/ (p)l(x-u) 1 . Further, 


f* Hi 

J. <*- 


Thus (see reference IB or 20), 


£>f (p)ip- a dp 

J. Cr-g) 1- 2tr J« (x-p) ' 


O (x, i) - O (x, a) where O !f,p) Is the Indeft- 
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which becomes, using the condition for Ci and relating x 3 and 
t) by the equation of ' the leading edge, Xi=2Motc— 1 = 
(2Wfc)-l, 


rfr 0 ) 


rCaVp 

P 


[* 


' 1 

2(3*r 0 


+ • 



( 68 ) 


By choosing the values of Oo, (ii, and bi given as follows: 


Mo 

Oo 

Ut 

bi=(a 0 /2)—(llp i )\ 


0.5 

56 

28 

26.67 ( 

(69) 

.8 

36 

27 

15.22 ; 



and placing the resulting expression for r(x s ) as given by 
equation (66) into equation (65), the values of — w/V 0 a 
shown in figure 31 were obtained. This figure demonstrates 
the accuracy to which the first two terms, of the series 
expansion for Iha,*) yields a constant value of w for the 
constants given in equation (69). 

The relation between circulation and lift has been derived 
and presented as equation (48). This expression can be 
written 

r d — r 0 

(=p 0 r 0 r+pocjp j r ‘^■v>(r 0 ,x 0 )dx 0 



where Xr,=x/c. In order to obtain a complete expression for 
the section lift, it is necessary to know the chordwise varia- 
tion of At p. Since equation (66) gives only the total vorticity 
and not its chordwise distribution, some assumption as to 
such distribution must be made. In lieu of this, the result 
presented in figure 27 suggests that for Iaige values of time 
the value of Aptro.Xq) used in the equation for section lift 




FicicnK 32.— Chordwise variation of circulation at end of early stage; A/i— Q.S. 


can be expressed by the product of A <p( °° ,x a ) and T(r 0 )/r( <»). 
In other words, for large values of r a , the shape of the chord- 
wise distribution of vorticity is the same as the two-dimen- 
sional, steady-state value. An indication of the accuracy of 
such an approximation is shown in figure 32 for a Mach 
number equal to 0.8 where the precise value of A^(r 0 ,ar 0 ) is 
compared with the approximation at t 0 = 4. 

Since T( “) = 5raT' 0 c/j8, the substitution of 

A $<1-0, X 0 )=A<p( a. , X 0 )r(ro)/T( co) 


gives for the section lift coefficient 

2r , 3 1 dT 
ei ~VoC + 2 IV dt' 


(TO) 


By means of equations (66) and (70) and the values. of the. 
constants given in (69), an expression for c tec can be written 
which is valid for large values of r 0 . This expression is 
somewhat cumbersome, however, and it is difficult to apply 
in subsequent analysis. Hence, it has been replaced by a 
simpler equation that is equivalent 17 to three decimal places 
for values of r 0 greater than about 10 for both Mach num- 
bers 0.5 and 0.8. These equations are given by the following: 
for A/o=0.5 

2t ["■< 4 1 44.218 1 , ' 

C ‘“ P L 3 5+2r„ (5+2^ w ' 


and for d/ o =0.8 18 


mi .in-uj. 

cbl A/i—OJS. 

Fpivai 31.— Variation of dewnwa^h at three-quarter-chord position. 


Ci 


2t r 

1 

1 P L 


1.736 


ll + 5ro/4 


70.83 1 

'( 11 + 5 ^ 4 )^ 


(72) 


ir in order to facilitate fairing the initial portion of the Indicial lift curre into the approximate solution obtained for the subsequent variation, the latter results were shifted slightly in terms 
of tV'/c. This did not affect the asymptotic behavior of the curve. 

* 1 * The value for c i a given by equation (72) disagrees slightly with the value given In the superseded TX 2403 (see Ibotnote I). The methods for calculating the results are somewhat 

different, and the value given here Is considered to be more accurate. 


2724S3— 54 2S 
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The equations for c m J, c t J, and for t q > 10 can be 
calculated from these expressions by means of the relations 
given in equations (64). Equations (71) and (72) will be 
considered valid for values of r 0 greater than 10 for both 
Mach numbers. Values of the indicial functions between 
r 0 =10 and the highest r 0 calculated by the exact method 
were approximated by interpolating between the two results. 
Tabular values of all the indicial functions calculated for 
the interval O<r o <10 are given in tables I and II and 
curves arc given in Part IV. 

TABLE I.— INDICIAL LIFT AND MOMENT FOR TWO- 
DIMENSIONAL SINKING WING. M„= 0.5 


I Vic 

fitlaf 2 * 

—2? Cm a '[* 

0.0 

1.103 

2 205 

.1 

.995 

L954 

.2 

.882 

1.632 

.3 

.772 

1.246 

.4 

.691 

.904 

.5 

. 650 

' .706 

.6 

.034 

.663 

. 7 

.618 

.002 

; .8 

. 604 

.561 

.9 

.588 

.540 

1.0 

.670 

.540 

LI 

.581 

.653 

1.2 

.007 

.572 

1.3 

.617 

. 000 

1.4 

.034 

.041 

1.6 

. 642 

.652 

1.0 

.669 

.067 

1.7 

.069 

. 082 

1.8 

.077 

.091 

1.9 

.680 

.699 

2.0 

.093 

.700 

2.6 

.724 

.727 

3.0 

.740 

.744 

8.6 ' 

.706 

.758 

4.0 

.781 

.772 ■ 

4.6 

.795 

.784 

0.0 

.800 

.795 

0.5 

.816 

.806 

6.0 

.826 

.810 

6.0 

.832 

.824 

7.0 

.840 

.832 

7.0 

.848 

.840 

8.0 

.854 

.850 

8.0 

.800 

-850 

9.0 

.806 

.804 

9.5 

.871 

.870 

10. 0 

.877 

.875 


TABLE II— INDICIAL LIFT AND MOMENT FOR TWO- 
DIMENSIONAL SINKING OR PITCHING WING. M 0 =O.8 


Wle 

to.*!** 

—Wem a 'jT 

^c,,'/3r 


0.0 

0.478 

0.956 

0.318 

—y • ; 

0.637 

.1 

.400 

. 920 . . 

.314 

.514 

.2 

.454 

.880 * 

.313 

.696 

.3 

.442 

.830 

• 21G 

- .581 

.4 

.430 

.750 

.326 

. 575 

.6 

.423 

.660 

-339 

.582 

.6 

.426 

.0Sf 

.366 

.696 

.7 

.433 

.632 

-307 . 

.808 

.8 

.442 

.620 ■ 

:a£5 

.621 

.9 

.451 

.615 

.397 

.633 

1.0 

.401 

- fill 

.411 

.644 

1.5 

.507 

.012 

.473 

. 6S8 

2.0 

.540 

.019 

.525 

.715 

2.5 . 

.581 

■ .628 

.560 

.729 

3.0 

.910 

.841 

. 698 

.737 

3.6 

.032 

.055 

. 026 

.746 

4,0 

.652 

.070 

. 648 

' .755 

E 4.0 

.670 

.084 

.085 

.765 

5.0 

.087 

.093 

. .681 

.773 

0.0 

.714 

.721 

.710 

.790 

7.0 

.738 

.743 

.730 

.800 

8.0 

.760 

.76! 

.700 

.821 

9.0 

.779 

.782 

.780 

.835 

10,0 

.798 

.798 

‘ .~79S 

.849 

00 

1. 000 

L000 

LOW) 

1.000 


INDICIAL FUNCTIONS FOR A TWO-DIMENSIONAL FLAT PLATE, Alt- 1.0 


The general results, obtained in the preceding section for 
the early stages of the motion and presented in appendix A, 
for the indicial loading over the sinking and pitching wing 
may be extended to the sonic case. Furthermore, the two 
intervals for which analytic results in a closed form were 
presented in appendix A now cover the complete time range 
since 0< r 0 < M 0 /(l-\-M 0 ) becomes 0<r 0 <0.5 and A/ 0 /(l + 
M 0 )< r 0 < Mo/(l— M q ) becomes 0.5<r o < “• Hence, by an 
appropriate limiting process, equations (A8), (A9), (AlO), 
and (All) become for 0<r a <0.5 > 


c, a =4 

Cm q '=-2+t 0 1 

C| /=2+r 0 l 


c « / = — (4/3) — (2/3) to*. 

and for 0.5 < r 0 < eo 

. • ■* ■ -- % 

c *«=W ir ) ^2V2r 0 — 1-f arc 


(73a) 


: COS 


To- 


ro 


l ) 


C *«'=— (2 M py^V2r 0 — 1+^1 -- y) arc cos 

<V = (2/ir) ^-V2 t 0 “1 + ( 1 + T | 2 ) arc cos T °~-~] 


r (73b) 


/=(4/3t) [ 14 - T g 3r % /2^1 + 

. : . o^) 


arc cos 


To- 


ro 


l ]J 


Since the magnitude of the functions in equations (73b) 
grows indefinitely with increasing time, the assumptions of 
linear theorj- are eventually violated. However, for moder- 
ate values of r 0 , these functions have the same order of 
magnitude as similar indicial curves for Mach numbers other 
than 1. These effects are illustrated in Part IV. 


INDICIAL FUNCTIONS FOB A TWO-DIMENSIONAL FLAT PLATE, ,U»-U 

The method of obtaining solutions for the indicial functions 
at supersonic Mach numbers parallels the development 
presented for the subsonic Mach numbers. The steady- 
state analog to the supersonic unsteady wing problem is a 
constant-chord wing tip with a supersonic trailing edge. 
(See figs. 33 (a) and 33 (b).) It is well known that the 
problem of finding the loading ovgr wing plan forms with 
all supersonic edges is one of the simplest in three-dimensional, 
lifting-surface theory. In fact, since the upper and lower 
surfaces .are non interne ting, the solution is determined by 
integrating sources within the Mach forecone. The analysis 
for c ta has already been carried out in reference 21. 

The analysis used to calculate the loading in terms of 
Xo=xfc and to=t/c over the three regions shown in figure 34 is 
outlined in appendix B. An example of the manner in which 
the loading varies with time over a sinking or pitching wing 
traveling at a Mach number equal to 1,2 is given in figure 35. 
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(b) 


(aj I Steady case, 
fb) Steady-state case. 

E-'iohre 33.— Boundary conditions for a two-dimensional unsteady wing moving 
at supersonic speed and the analogous three-dimensional, steady-state wing. 

The expressions for the indicial lift and pitching-moment 
coefficients are given analytically in appendix B, and plotted 
in Part IV. 1 1 can be shown that the results given in appen- 
dix B reduce to the expressions given by equations (73) 
when Mo is allowed to approach 1, so that there is no dis- 
continuity in the theory in passing through the sonic range. 

SUPERSONIC STEADY-STATE UPT 

Since the next section contains a problem involving a 
complicated acoustic plan form, it may be helpful to consider 
first a problem involving a very simple acoustic plan form 
but otherwise similar to the subsequent analysis. Hence let 
us inspect, using Kirchhoff’s formula, the problem of finding 
the steady-state loading on a two-dimensional flat plate 
traveling at a constant speed (equation (15)) . 

Since the upper and lower surfaces are noninteracting, 
we can use the special form of Kirchhoff’s formula given as 



Fitir aK Sf.— Regions used In analysis of supersonic unsteady wins. 


equation (28). In the plane of the wing this equation be- 
comes 

1 f f Wudxidyi 

? 2xjs. J V(z — *i)*+(y— 2/0* ^ 

Figure 36 shows the positions of the wing in the xy plane. 
The wing has constant speed for t>0. The point P(x,y) 
is chosen on the wing and ahead of the wave which started 
at time zero; therefore, P{x,y) lies in the region which has 
attained its steady-state value. Further, the value of vr u 
is constant over the acoustic plan form in such a region. 
This constancy reduces the problem to one of integrating 
[Cz— Ti)*-f-(y— yiY\~ ul over the ellipse representing the 
acoustic plan form. 

The equation for S a can be determined from equations 
(20) and (21). In this case, equation (21) becomes simply 

Xi= —MqT 

Eliminate r between tills equation and equation (20) and 
there results in the 2=0 plane 


(x-x 1 )*+(y-2Ii) i =(i+^y 


That equation (75) is the equation ot an ellipse with one focal 
point at a :,y can be readily verified. It is more convenient, 
however, to change to a polar coordinate system with origin 

at P. Hence set 

x—Xi =r cos 8 

y—yi=rsm8 

dy l dii=rdrd8 

Then equation (75) becomes 

J/,>r= J/ 0 t+ar— r cos 8 


x -pMj, 
J/o+eos d 
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(a) Sinking wing. 

(b) Wing pitching about Its leading edge. 

Figure 35.— Variation of two-dimensional Indidal load distribution with percent 
chord for a Mach number equal to 1.2, 



Figure 3i.— Xmustlc plan form for point In steady-stato region of two-dimensional un- 
steady wing flying at a supersonic speed. 

and, therefore, equation (74) becomes simply 

i it - 

•r ' x± aw 

W u C 2 * Jn rMt+CO*t 

■ , 2i]» H dr 

'' r 2r de 

- 2 ^(2+Mo0j o Mo+ - cos g 

The integral is not difficult to evaluate so 

and, finally, by equation (12) 

Ap 4 / w u M 0 \ 4 a , 

S FoA/a'v -y/MT^l) y r W~l 

which is the familiar Ackeret value for the loading on a 
two-dimensional flat plate. The lift coefficient,, of course, 
follows immediately as 


c,= 


-yW-1 


(78) 


INDICI AL LOADING FOE SINKING TBIANGULAR WING WITH SUPERSONIC 

EDGES 

The geven regions. — The analytic expression for the indi- 
cial loading over the triangular wing has a different form in 
each of seven regions. These regions are determined by 
the positions of tkj various wave fronts relative to the wing 
plan form" (fig. 37). For t<C . 0 the wing is motionless, its 
leading edge lying along hues represented by the dashed 
lines in’ figure 37. The trailing edge is considered to be 
supersonic and hence its position is immaterial, since the 
solution can ho cut off wherever desired. At i=0 the wing 
starts suddenly to move, and for 0, travels forward at a 
constant speed T 7 0 . After a certain time t has elapsed, the 
wing has traveled to a new. position, also shown in the figure. 
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In this same interval of time, pressure impulses have trav- 
eled out in spherical waves from every point of the region, 
which the wing has occupied. The trace on the wing of the 
sphere starting from the wing apex at 1=0 forms the exter- 
nal boundary of region 7. The area outside this circle and 
within the traces of the cylindrical waves (the envelopes of 
the spherical waves) generated by the leading edges at f=0 
forms region 4. Region 5 is formed by the overlappping of 
these cylindrical waves, and the solution for loading within 
it can be found by a suitable superposition of the solutions 
for regions 3 and 4. Region 1 lies between the cylinder 
trace on the wing and the leading-edge position at time f; 
the loading in this region cannot be affected by the manner 
in which the wing started its motion since it lies outside the 
starting cylindrical waves. Hence, the loading in region 1 
is the same as that on a swept wing flying at a steady super- 
sonic speed. The solution in region 2 can also be obtained 
from steady-state lifting-surface theory, hut, whereas in re- 
gion 1 the field is two-dimensional (i. e., invariant with dis- 
tance measured parallel to the leading edge), in region 2 
the field is conical. Region 6 is formed by the overlapping 
of regions 2 and 4. Finally, region 3 is that area com- 
pletely unaffected by waves from the wing edges. In the 
following subdivisions the analysis of each of the separate 
regions will be discussed. 

Region 1: The loading in region 1 of figure 37 is equal to 
the loading on a two-dimensional flat plate moving at a 
constant velocity given by the component of stream velocity 
normal to the leading edge of the triangular wing. Since 
this component is supersonic, the loading is of the Ackeret 
type and is given by 

/ Ag\ 4a„ 

'?»/i A. 

But since F*=F, cos A where A is the angle of sweep (see 
fig- 38), 

qn—qo COS ' 2 A 

a m = a sec A 
d/,=il/o cos A 
j9,=yil/ 0 * cos 2 A— I 

and 

/Ap\ __ 4 a 

V?o/i V-Mo s ~sec 2 A 




Figure 38.— Equations of lines used in. the analysis of the triangular wing with super- 
sonic leading ed ges . 


Finally, if ctn A =m 



where jS = V-Vo 2 — 1 . 

Region 2: The steady-state loading on a triangular wing 
with supersonic edges has been given by several authors (see, 
for convenience, reference 22) so the expression for the load- 
ing in region 2 can be written immediately for the coordinate 
system shown in figure 38 as 




sr+arc sin 


fihny— (j+A/pt) 
P[m{x+Mat)— y] 


« m P 2 my+(x+Mot) ) 
ftmix+M^+y]] 


(80) 


Region 3 : Since region 3 is unaffected by the edges of the 
wing the solution for the loading therein can be written as 
in reference 21 



Region 4: The solution for loading in region 4 can be 
obtained from consideration of a two-dimensional wing 
starting from rest and moving with velocity F, normal to 
its leading edge. This problem has been treated in refer- 
ence 21 and the solution written there can be written for the 
right-hand side of figure 37 as . . . 



M m x,+t V-AT » 2 1 

x.+Af. r 



where the notation, as defined by figure 39, is 


Fiuc*e 3“.— The seven, regions used In the analysis ot the triangular wing with supersonic 

trading edges. 


x x =x cos A — y sin A 
y K =x sin A-j-y cos A 
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and since 


cln A =m, sin A= 


V 1 +m 2 


Then 


cos A= 


m 


V 1 


x * : IZ 


Vl+m 2 

x+my 

y » — — _ 


/ 1 -fm J 


The equation for ioadhig now becomes, in the coordinate 
system of figure 38, 


/Ap\ 4 am 

\ ?o )t ir-y /jS*m s — 1 


arc cos 


wMofma;— |y[H- <(1 4-wt 8 ) 
VT+m^Cmr— \y\ +m.W^) 


V)3 2 m s — 1 /V , . mi— 

— : — I - 7 + arc sm — p= — - 

0 V 2 tVl+m* 


mil/, 


(82) 
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Figure 40. — Solutions superimposed to obtain equation (83). 


solutions for regions 2 and 4 (equations (80) and (82)) and 
subtract the solution for region 1 (equation (79)). There 
results 


\ 9o/ 8 ir- 


4am 


Vl 8 ! m>- ll 


arc sm 


P*my — (z + MbQ 
/3[m(z + Mo<)— j/l 


. jS s my + (:r + MV) , 

areS 1 %- [m("+J^+ri + 

arc cos ?^«!Hd3d)+Al±22- 

(m a- — 1 3 / 1 + m Af d) Vl+m* 


Vj3 2 m 2 — 1 /tt in ma- iKTYl 

milfo ^ 2 + iVf+m* 


(84) 


Region 5: Tlie solution for loading in region 5 can be 
obtained by superposition of the solutions for regions 3 and 
4 . (See fig. (40).) If the solutions for the two sides of 
region 4 (obtained from equation (82)) are added, the result 
gives twice the required value of w u on the wing, as well as 
undesirable pressures off the wing. However, subtraction 
from this sum of the solution for region 3 (equation (81)) 
reduces the downwash w u to the proper value, and also 
cancels the excess pressures. The resulting expression can 
be written 


Region 7 : The solution for the loading in region 7 can be 
obtained by means of equation. (28). The analysis used in 
finding the solution in this region is not difficult but the. 
algebra, is rather involved. It is useful at this point to intro- 
duce polar coordinates (see fig. 41) such that 

x—Xi=r cos fl'] 


y—yi=r sin 8 f 


(85) 


<hidy L =rdrd.B J 


/Ap\ 4am r cog TOM 0 (mz-y)+(l-tm*)t ^ 
\9o/ 6 irV^m*— 1 L Vl +m 3 (mx—y+mM j) 

arcco, 7 M « ( " I+y) V 1 + ’f | ‘+ 

V 1 +m\mx+y+mAIJ) 


1 ( arc sin 

mM 0 V 


mx-ty 

t-Jl+m* 


arc sin 


mx—y V 

t-y/l +fflV. 


(83) 


Region 6 : The loading in region 6 can also be calculated 
by superposition. To find the loading in this case, add the 


From oquation (85), equation (28) can be written in the 
form 

•“Tf/J** (8C) 

The acoustic plan form for points in region 7 is the region 
bounded by three curves as indicated in figure 42. The arc 
between fli and 9 t is determined by eliminating T between the 
equations 

r*=(f— T) s 
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(the equation for the inverse sound waves) and the equation Using these expressions) equation (86) reduces to 
for the left leading edge 


2/i = — m(2,+Af<,T) 

The arc between 0* and 8 { is found by determining the acoustic 
intersection (eliminating T) of the right leading edge with 
the inverse sound wave; and, finally, the arc between 0* and 8 2 
is given by the equation r=t. The equations of these arcs 
can be written in polar coordinates as 

r= — e l <e<e i 
m cos fl-hm-Uo+sni 6 ~ 

r=(;0 a <0<0, 

m(x+A(g£)—y f 0 < <9 < t 

m cos 0-j-mM o — sin 6’ (0i<0<2x 



Fiqctjh 41.— The rt coordinate system. 


Wing at t~o 

Wing now 

B\ -8, -Zn 


_ Voa f* 1 m(x+M^)+y , , 

V 2x Jty m cos 0+m2lio-f-sin0 a 

Vo* P* 1 m(x+AU)—y dfj 

2r J» 3 mcos0+mM o — sin0 

and taking the partial derivative with respect to f (to deter 
mine the loading according to equation (12)) one finds 18 


/Ap\ 2am f** d& 

\ 2o/t~ it J*! mM 0 +m cos 0+ sin 8 

2 am r'x dd 

x mM 0 +m cos 8 — sin 8 


2a f* i 

'tmJh 1 


In evaluating this equation, the following integral is used: 
for — x<0<x * 

C de 2 

J mAf 0 +m cos 8± sin0 1 

m(M»-l)ta-(»/2)±l 

Vi s’m 1 - 1 


Since equation (8S) is valid only in the interval— x<0<x, 
care must be exercised in applying it because the angle 0, 
may be greater than x (as in fig. 42). In case x<0„ it is 
convenient to introduce the angle 0, , =0 J — 2 t. The expres- 
sion for Ap/g, can then be written in two forms, according 
to whether 0, is less than or greater than x: 


for y> 0, 0,<x 


/ Ap\ 4am T . , m(4/ 0 — 1) tai 

V ?o/7 xV^*m s — 1 [_ V^m 1 

m (M 0 — 1) tan (8J 2) -f-1 , ‘ 
arc tan — - — - — , — — + 


— 1) tan (0^2)— 1 


arc tan 


arc tan 


V/Sm’—i 

m(l/ 0 — 1) tan (0^2) + 1 
V^’m*— 1 

m(M a — 1) tan (0^/2) — 1 ~ 
Vj3*m* — 1 




for y> 0, x<0, 


/Ap\ 4am r fBn m(xl/ 0 — 1) tan(0i/2)— 1 

V?o/7 xV(9 J m J — 1 L V/3*m* — 1 

m (J/ 0 — 1 ) tan (0^2) -f 1 

arc tan — - — , == = —*■ h 

Vj3*m 1 — 1 


arc tan 


m(ilf 0 — 1) tan (0a/2)-j-l 

V/?*m* — 1 


tan m(1I ° ^) tail( y /2)— 1~| + 2a ( ^_ fli+2jr) 

Vj3*m*— 1 J T Ai 0 

(89b) 


Fig urr 43.— Acoustic plan form for point In regf on 7. . | 

!, Tb« limits fft, ffi, and tf* are all functions of t but In moving the partial derivative through the Integral sign the terras Involofng 'Vhftt, and all cancel one another. 
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where 


= arc cos 
(0 «?!<*) 
(? 2 =arc cos 

(O<0J<7r) 


— Moyt+jx+Mfti-Jix+Myt ) 2 — ffV 
y 2 +(x+Mj.y- 


m(y+mx) — ^{l+m 2 ) tf—jy+mx)* 
(1 +m 1 )t 


—fn(y—mx)—^{l+m 1 )t i —(y — mx) i 

0,— arc cos (i +m*)f 

6) 3 7 = (9, — 2ir (if 7T < <?a) 


The limitation on 6 S can be given both an analytic and 
geometric interpretation. Thus, equation (89a) applies for 
0 <m(x+t)<y and equation (89b) applies for 0 <y<m 
(as+t). These regions are shown in figure 43. Because of 
the geometrical symmetry about the x axis, equations (89) 
suffice for the determination of loading throughout region 7. 

An isometric drawing of the load distribution on the right 
panel of a triangular wing with supersonic edges is shown in 
figure 44. The positions of the spanwise sections were 
chosen so that each of the regions 1 through 7 is represented. 
It is to be noted that the results for region 7 show no unusual 
characteristics and, in general, the distribution is similar to 
the steady-state loading on a triangular wing. 



INDICIAL LIFT AND PITCHING MOMENT ON A TRIANGULAR WING WITH 
SUPERSONIC EDGES 

The indicial lift aud moment could be obtained, of course, 
by integrating load distributions calculated by the method 
presented above. However, it is far simpler to use the 
methods outlined in the section entitled “Boundary-Value 



Fionas 41.— Distribution of IndlcM loading on right panel of suporsonlc-odged, tri- 
angular wing. . 


Problems Involving Noninteracting Surfaces” in Part II. 
In particular the second of the two methods outlined therein 
will be applied. 

It was seen in the above-mentioned section that the lift 
and moment coefficients for a supersonic-edged triangular 
wing could bo found by solving a related steady-state lifting- 
surface problem. For the cases of indicial sinking and pitch- 
ing, the boundary conditions for this related problem can 
be found readily from equation (34). Since y r and y, 
become, for s=0, the right and left leading edges, respectively 
the boundary conditions are given by 

Gf), a = ~ 2V o am ( xJrM ° t ') (90) 

for the sinking wing, and 

(M)..r ~ 2kn{x+m)2 C91) 

for the pitching wing. In addition, it was shown that the 
lift and moment coefficients for the wing in unsteady motion 
are obtainable by integration of the quantity For 

convenience, this quantity is expressed as a chordwisc load- 
ing factor in terms of the following notation: 


ss.i i 
II 

& A© 

— dy (sinking wing) 
-i Q o 

(92a) 

mM 0 V 0 

1 —dy (pitching wing) 

(92b) 

.2 M 0 
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where 

It is found that (equation (35)) 

(Ji-KWi) dx v dt\ 


1 " — El 


Tit c ytj.j 
4 m* d 


= 4 r C jx^M^Ydx.dt, 

1 tts “ cM.J J yff—t,)*— (x— ard* 


(93) 

(94) 

(95) 


The lift and moment of s inking and pitching wings will now 
be obtained by the use of P 0 and Pi. 

Wing with constant angle of attack (si nk i n g wing). — The 
solution for the chord loading on a flat, supersonic-edged, 
triangular wing starting from rest at f=0 and flying at a 
constant speed and angle of attack is given by equation (94). 
With the transformations 


X—Xi=Xi 
t — t\ = tt 


this becomes 


f S * j— 

V-1 dU<rJ \tY-X-r 


This integral can be evaluated and gives for x <£ (region A) 

P 0 =4 (96a) 

for —t<x<t (region B) 

t + AI 0 x J 


, _4 r yV-x* 

° IT _X -f- 4/ tji 


-arc cos 


-i/of 

and for r<— f (region C) 


B) 


-1/o 

P 


+ ^ arccos ^J (96b) 


it—g- 


(96c) 


where the regions are shown in figure 13. As has been 
pointed out, equations (96) can also be obtained by inte- 
grating the equations for the loading given in the preceding 
section. These integrations were carried out (in some regions 
numerically) and the results were found to agree with those 
of the present analysis. 

It is now possible to write the indicial functions C La and 

C m ' in the form 

« 

O — A/q/ 

= <?V m(x+M<f)P a dx (97) 

aj n — M 0 t 

r* w # " mix+M'tyPodx (98) 

OCn-Mrt J —Mat 

where c 0 is the root chord, S is the wing area (equal to mc a l ) 
and the prime indicates that the pitching moment is measured 
about the apex, the positive moment being one which causes 
the trailing edge to sink relative to the apex. 


Interval 



Figvre 45.— Time Intervals used fa exprrsfa* lift and moment on unsteady, super- 
sonfc-edged, trijnjuHr wing. 


Combining equations (96) and (97) one finds for the first 
interval shown in figure 45 

c ‘-jfeiJ-l i f 0(jr+JWl+ 

; cos (-f)+ 


^ arc cos dx+ P* ^ 4(x+SU)dx \ 

P -x+MJJ Jt > 


This equation integrates to give, if fo=f/c 0 , 

for 0<f tt < ,/ lT ( first interval) 
xUo + 1 

(99a) 

Similarlv for v-A-r^^TT^T (second interval) 

MQ-rl A*u l 




arc cos 


A/(/o 1 


— arc cos {M (i —t 0 p x )+ 

2 3 ~y° - VtfB 

TTAZq 

and for . - 1 — -<t 0 (third interval) 
At o A 

r -i 

Ll ~p 


(99b) 


(99c) 


In the same manner the values for C m J in the various 
intervals can be determined by combining equations (96) 
and (98). There results 

for 0<f a < , r\. (first interval) 

/vAoi" l 

c -.'~m( 1+1 P w ) 


n r = 

L "a 


(100a) 
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(second interval) 


(8 -M„<o-M 0 V- 2 f 0 2 )V«o 2 - (1 -MW+ ' 

i (2+Mdh?) arc cos — -+^ ft rc cos (Af 0 — iS^J 

(100b) 

and for —* a ( fc hird interval) 



(100c) 


Numerical results will be presented in Part IV. 

Wing with linear angle -of-attack variation, (pitching 
wing). — The solution for the average load on a flat, super- 
sonic-edged, triangular wing flying at a constant speed and 
pitching at a uniform rate 6 about its apex is given by 
equation (95). With the transformation 


X—Xi—Xi 


this equation becomes 


P i= 


4m 2 d 


jrs 2 dt 


SJ 


( m Xl M °* 2 ) 


■y/tf-xS 


— dttdx % 


and the evaluation of this gives (for the intervals defined in 
fig. 13) 

for x< t (region A) 


for — t<x<t (region B) 

„ in t’+ix+Mjy ( x\ 

p ’—*h (i W '*"”!-!) 

Ala t+Jl/oJ . 1 3 2 -{-2 AdJ- 

T“ eooa ^+m + 2 -&+m 

and for x< — t (region C) 


(101a) 


+ 

Vf^-aT* 




4M« 


(101b) 


(101c) 


The equations for the indicial functions Cl * and C Mq can 
be obtained from the equations, 

f£r°<-s£n. f-L? <>&+»wp,dx d02) 

MV 

/FT (wo) 

MV 


where the primes indicate the wing is pitching about and the 
moments are measured about the leading edge. 


A combination of equations (101) and (102) gives for the 
lift coefficient 

for 0 <t n < ,> y -r- (first interval) 

/■ nio+l 

(104a > 

for *, ■ r- <t< tt — r (second interval) 

~M 0 —l 1 

°-:~m [(I w - j **’ +s) arc cos - t -i+ 

^ arc cos (M„—pHd+ 

(|<o ! +| A Mofo+y|)v'fo s - (1 “ W] (104b) 

and for 3/Trj— ^ interval) 

C h '= 8/3j8 (104c) 

Similarly a combination of equations (101) and (103) yields, 

for the pitching moment about the apex, the results 

for 0 ry (first interval) 

0*/=-^ [l +io*-| V(1 +43/0 2 )] (105a) 

for -^.F— . <t a < ^ (second interval) 

[[ 2(1 +«-? w + 4W ] “ rc 

^arccos ( M 0 -fPt„)+ 

P 

[42-22 AU « + (2 AI a 2 + 3)£ a 2 + (2 Af 0 ! + 1 3)il/</o 3 ] 
V't 0 2 -(1 --W 0) 2 j ( 1051> ) 

and for -j j y < t 0 (third interval) 

<V=-f (105c) 

Numerical results are presented in Part IV. 

INDICIAL LOADING FOR SINKING TRIANGULAR WING WITH SUBSONIC 
LEADING EDGES 

The six regions. — As in the study of supersonic-edgcd 
triangular wings, there are also in the case of triangular wings 
with subsonic leading edges various regions in which the ana- 
lytical form of the loading equation is different. Figure 40 
shows the regions into which the subsonic-edgcd triangular 
wing can be most conveniently divided, where the trailing 
edges are again assumed to be supersonic and the solutions 
are cut off appropriately. Most of these regions have 
counterparts on the supersonic-edgcd wing shown in figure 37. 
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To begin, with, region 6 lies within the spherical wave which 
startecT&t f=0 from the wing apex. Region 1 is within the 
cylindrical wave which was started at f=0 by one edge of the 
wing, but outside the wave started by the other edge. Re- 
gional is the area formed by the overlapping of the two 
cylindrical waves from the opposite edges, but outside the 
region influenced by the reflection of one of these waves on 
the opposite edge (secondary wave fronts shown in the 
figure). Region 5 is the area between regions 4 and 6 where 
the flow is influenced by secondary (and higher order) wave 
reflections. Finally, regions 2 and 3 are similar to regions 2 
and 3 in the super3onic-edged case; region 2 being that unin- 
fluenced by. the starting phenomena and therefore having a 
loading already at its steady-state value, and region 3 being 
that which is unaffected by the disturbances emanating from 
the edges. 



Fftirai W.— The six regions used tn the analysis of the triangular wing with subsonic 

leading edges. 


Region 1 : The solution for the load distribution in region 1 
is the same as that for a two-dimensional wing starting 
suddenly from rest and moving with a steady subsonic 
velocity U* normal to the leading edge. A solution to the 
latter problem for the initial part of the motion is presented 
as equation (A6(b)) in appendix A. In terms of the normal 
components of velocity and distance, therefore, the loading 
coefficient for the right-hand side of figure 46 can be written 
immediately: 


(*p\ = 

8u? 0 ( A/* / 


- TV m Zin\l+M u \ 


t—x n 


— [-arc tan 

ft 



The equations which relate the normal components to 
those in the free-stream direction have already been given 
in the section on region 4 of the supersonie-edged wing. 
Use of these relations leads to the following expression for 
loading in region 1 (in the coordinate system of fig. 47) : 

/Ap\ _ 8a ( mAI Q lt^\ + m,*+\y\—mx . 

Wo A VmA/o+vT+m* V mM 0 t-\y\+mx _l ~ 

/ mM a t—\y\+mx \ 

arc tan ■■■ ■- ■- , 4 (106) 



Figs:** 47.— Equations of lines used fax the analysis of the ftfcwgulsr win* with 
subsonic Ieadln* edges. 


Region 2 : The loading on region 2, being the steady-state 
loading on a triangular wing with subsonic edges, is well 
known. The solution for region 2 of figure 46 is therefore 
given by (see, for convenience, reference 22) 


A p\ 4 gm z (ir+M)t) 

. 2 o / s E\'m 1 (x + M a £)*— 


(107) 


where E is the complete elliptic integral of the second kind 
with modulus Vl — 0 2 m*. 

Region 8 : The loading in region 3 follows from reference 21 
and is 


Ap\ 4a 

. 2o/j Af 0 


(108) 


Region 4: The loading in region 4 of figure 46 is calculated 
by superposition, just as the solution for region 5 of the 
wing with supersonic edges was obtained. The solution 
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in region 4 is the sum of the solutions for the right and left 
halves of region 1, minus the result for region 3. Thus 


/A;A _ 8a mM 0 I /tyl +m*— y— mx . 

\ go A irA/ 0 LmA/ 0 +Vl +m* W ™.M 0 t+y+mx 


_y/fVl +m 2 +y — mx 


mMot—y+mx 


arc tan 




j + arc 


mM 0 t— y+mx 
/l+m 2 +y— mz 


mMpt+y+mx 

/l+m 2 —y—mx^~ 

(109) 


Begions 5 and 6 : In these regions the exact solution for the 
loading has not been determined. As was shown in the 
section on homogeneous boundary-value problems, such 
solutions would require the solution of a three-dimensional 
elliptic-type partial differential equation. A later section 
will contain an approximate solution for these regions when 
the wing plan form is slender. 

Discussion. — An isometric drawing of the load distribution, 
for the regions in which it is known, is shown in figure 48. 
Comparing the results for the loading on this wing to the one 
with supersonic edges (fig. 44), it is apparent that the prin- 
cipal difference in the two distributions is in the behavior 
around the leading edges; the loading being finite at the 
supersonic edge, whereas it becomes infinite at the subsonic 
edge. In view of the known steady-state results lliis differ- 
ence was to be expected. Elsewhere the loadings are quite 
similar. 


Traces of regions shown m 



Figure 43.— Distribution of indlcfal loading on right panel of subsonlo-edged tri- 
angular wing. 


The results presented in equations (IPG) through (109) will 
next be examined in a different light. Choose a given span- 
wise sqctjqn on the w'ing and watch this section as time pro- 
gresses; from £= 0. This amounts to fixing the axis on the 
body and can be accomplished simply by using the quantity 
8 introduced in equation (93), 

8=m(x+Mot). 

It is clear that « is the semispan of a given spanwiso sections 
and that if equations (10G) through (109) are written in term, 
of s, y, and t, for a fixed s they represent the variation of 
loading on a given section as time progresses. 

If the notation is further simplified by introducing the 
parameter P e where 


&•— 


1 

jmA/o+Vi + w J 


- ( 110 ) 


equations (106) through (109) can ho written in the following 
way: 


/A p\ 4 am 8 

\2oA E\s 2 —y 2 

f&p\ 4g 

\ 2 o / a M, o 


( 111 ) 

( 112 ) 

/ 

(113) 




— 8 


+ arc tan 


arc 


«+y 

tan V(£/&)-y-e“ 2) 


/ "V-y ' 


■+ 


(114) 


The load distribution across any section is given by equa- 
tions (111), (113), and (114) from the time <// 8,=0 to (£//?,) 1 , 
where the term is equal to 2s or s/m(Mp+ l)P„ which- 
ever is smaller. (At t/p e —2s the secondary waves shown in 
fig. 46 have just reached the spanwise section, and at t/p t — 
s/mjS e (Mo+l) the spherical wave which started from the 
apex has just reached the spamvise section.) From {tjp e ) , 
to (f/j3 e ) a =s/ro&(A/ 0 — 1), the loading has not been determined 
and from f/&=(f/i8«) a to t = co the loading is the steady -state 
value given by equation (112). Figure 49 show's this initial 
and final load variation plotted as a function of the parameter 
tfp t . At the beginning of the motion the loading is constant 
across the span, but this type of distribution is quickly 
modified arid the shape of the curve tends toward the steady- 
state loading given by equation (112) and shown in the figure 
as the distribution at £//3«=(£/0«)s. In fact, if the span is 
crossed first by the secondary w’aves rather than the spherical 
wave, when this span has traveled a distance such that 
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f//9 4 =2s, the expression for the loading given by equation 
(1141 becomes 

( \p\ _lQm(i,8a 


( ■±p\ _ lbmp,8i 
9oA_ w V**— y 


which differs from the value given by equation (1 12) only by 
a constant of proportionality. Both before and after the 
time f//3«=2a the shape of the loading curve varies from the 
simple type represented by equation (115), but the trend is 
established. 

The average chord loading factor P 0 , introduced by equa- 
tion (92a). can now be determined for certain regions. 
Hence, if the notation 

t—Us \ 


p 0 =^- # - f* — d y 


is adopted, there results for the early part of the motion, 
that is, for 0<r//3,< (t/| 3,): 

/’.>= 2 ( 2 -£) + 4m^.4/ 0 W(S.) (117) 

Equation (117) was derived by integrating equations (111), 
(1 13), and (1 14) . For values of r//3, > I/jS,m(M 0 — 1)=( ~ ) 

.Ptf / 2 

equation (112) is valid. Hence for-£->( 

Pt \P«/i 


»— e 


(118) 


Figure 50 indicates the magnitude of this average load for 
both large and small values of r/0«. Notice that for small 
values of r//3« it is sufficient for the establishment of the 
curve to specif 3 r the parameter but for large values 

an additional parameter must be given (such as M a in the 
figure). Notice, further, that in spite of the large variation 
in the distribution of the loading, as shown in the previous 
sketch, the average value P a varies linearly throughout the 
intervals considered. This result is s imil ar to the one 
obtained for triangular wings with supersonic edges. 



0 4 8 12 18 20 

r/P. 


Figcsk 50.— Variation of P t with t % at a fixed spanwis« section. 

INDICIAL LOADING ON VERY SLENDER TRIANGULAR WINGS 

In Part II of this report it was pointed out that if the wing 
is slender (i. e., has a small ratio of span to chordwise length) 
the governing partial differential equation simplified to the 
form 

Vtt — <Pn~~ (119) 

where the Independent variables refer to a coordinate system 
that is either fixed on the wing or is fixed with respect to the 
still air at infinity. The boundary conditions that apply 
when the axes are fixed on the wing will now' be considered 
in some detail. (In Part II the variables x t , y u z x and t t 
were used -to denote this coordinate system. In order to 
avoid a cumbersome notation, however, the subscript 1 will 
be deleted in the following.) 

Just as in the previous sections of this report, consider a 
triangular wing which is at rest for f<0, starts suddenly to 
move at a forward velocity equal to H 0 at f=0, and continues 
at this same velocity for C>0. In this case, the triangular 
wing can be considered as having a finite chord, since span- 
wise sections act independently in slender wing theory, and 
the final integration in the x direction can be stopped at any 
desired chord length. It should be emphasized that in this 
case, Ho may be either subsonic or supersonic. A section in 
the spanwise direction, as for instance section A A in figure 
51 has a trace in the yt plane which is a narrow rectangular 
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Kiouus JL— Slender win* In zy plane and trace of section In yl plane. 


strip along the t axis. Since equation (119) has been derived 
on the assumption that the velocity gradients in the y, z, 
and t directions are independent of the gradient in the x 
direction, the boundary conditions along the strip shown in 
the figure are independent of those on other strips corre- 
sponding to traces from spamvise sections along the wing. 
Hence, the problem is to find a solution to equation (119) 
which will make <p, constant over the strip (since a given 
spamvise section experiences constant downwash whether the 
slender triangular wing is undergoing sinking or pitching 
motion) and at the same time will satisfy the other condi- 
tions listed under equation (2). In the lifting-surface 
analog this corresponds to the problem of finding the velocity 
potential over a flat rectangular wing of low' aspect ratio 
situated in a free stream moving at a Mach number equal to 
sj2. Solutions to the latter problem can be obtained by 
various techniques, and so the procedure will be first, to find 
the potential for the steady-state, flat, rectangular wing, atid 
then, by analogy, to convert tills to the solution for either the 
sinking or pitching slender triangular wing. 

The steady-state, lifting-surface problem. — Lifting-surface 
solutions for the loading on a rectangular wing traveling at 
supersonic speeds have been developed for regions 1, 2, and 
3 of figure 52 (by Busemann and others), and by means of 
these solutions the load distribution on a spanwise section of 
the triangular wing can be determined to a time necessary 
for sound to travel that span length. For however, 

the solution becomes considerably complicated by the in- 
creasing number of reflections from the edges. Reference 24 
gives solutions for the loading on a rectangular wing in 
region 4 and indicates methods for extending the solution to 



x 


.Fiouke 52 — Regions used In the discussion of the low-aspeot-ratlo rectangular wing. 

regions farther along the wing. Already in region 4, how- 
ever, the expression is cumbersome and in higher-numbered 
regions the expressions become difficult to manipulate. 
These methods, therefore, will be discarded in favor of a 
more approximate hut simpler analysis. 

If x is the distance along the chord, y the distance along 
the span, and s the semispan, then the solution for regions 
1, 2, and 3 of figure 52 can be wTitten (for convenience, see 
reference 22): 


Region 1 


Region 2 


A W 

Ho 


Ap 8w 

— ==- arc tan 

2o t Vo 


Vr-s-Hj 


Region 3 

— = - — - ^ fare tan -*/ — + a rc tan -*/— - — ~~ —'j 

2o *-To\ \ x—e—y ' V x-s+y 2/ 


(120c) 
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As x increases (i. e., for higher-numbered regions in fig. 52) 
it is reasonable to assume that the spanwise variation of 
loading is virtual^' invariant with x, that it is “smooth," 
and that it falls to zero at the side edges. Assume, there- 
fore, that the loading is given by the relation 


<io 





( 121 ) 


Spanwise, this has the variation shown in figure 53; chord- 
wise it is as yet arbitrary. To fix the chordwise distribution 
the value of f(x/s) will be determined so that the vertical 
induced velocity along the center line is constant and 
equal to w. 

The solution to this somewhat artificial problem approaches 
the exact solution to the steady-state lifting-surface problem 
for a flat rectangular wing along sections far behind the lead- 
ing edge; closer to the leading edge it only approximates the 
exact solution; and, of course, in the vicinity of the leading 
edge it will be least representative. But-, on the other hand, 
the exact solution is known in the vicinity of the leading 
edge and it turns out that the solution of the problem posed 
above forms a reasonable continuation over the remainder 
of the wing. 

The velocity potential for the problem which has been set 
can be readily expressed in terms of an integration of ele- 
mentary horseshoe vortices over the plan form. Since the 
Mach number equals -^2, then according to reference 19, 


v r 0 ~ |’ f ( x — x l )(&plqg)dx l dy l 

* -It JaJ Xi)*— Vi*— z 1 




Fiours 53. — Spanwffe loading and region of Integration used in th* analysts of the low-espect- 
railo rectangular wing. 


where A is the area on the wing within the forecone from the 
point P(x,y,z), at which <p is to be determined (the shaded 
area in fig. 53). 

The simplification of the last expression is given in refer- 
ence 19. The result is the integral equation 

1 =/(i?)+— f /(>?i)G r (’j — (122) 

where ?=x/« and G is given by 

l <n~vi 


U — 77t < 1 


The modulus of E\ is k\ and the modulus of K t and Et is k t . 

The solution of equation ( 122 ) for f(jj) is not difficult 
when numerical methods are used. For intervals of 171 equal 
to 0 . 2 , the result is given in tabular form in table III, and 
also in figure 54. As mentioned previously, the function 
f{rj) determined by this approximate method will be least 
representative of the exact solution in the region near the 


G(.v—vi)=-{ 


E l 

Er-A-k&Kt 

£2 


kr- 


1 


V~Vi 


; k i =ri — ij l 


leading edge. However, since the exact solution is known 
for x/s=? 7 < 2 , a compromise can be effected. If the span- 
wise average of the loading is calculated from equations 
( 120 ), it is’ found to be a linear function of 17 , starting at 
4/y for 17=0 and falling to zero at 17 = 2 . On the other hand, 
the spanwise average of loading given by equation ( 121 ) is 

— 7 r-fr/(j 7 ). Therefore, an improved solution for /( 17 ) consists 
v 0 

in taking/(rj) as the value given by 


/0l) = ~(2 11 ) 

in the region 0 <t?< 2 , and then fairing this curve into 
that given by the solution of equation ( 122 ) for i 7 > 2 . The 
two curves are shown in figure 54. 

By using the results listed in table III, the loading over a 
low-aspect-ratio rectangular wing flying at a Mach number 
equal to y 2 can be estimated. Of particular interest is the 
damped oscillatory nature pf.the load, as shown in figure 54, 
falling to zero at one span length behind the leading edge 
and taking alternately negative and positive values beyond 
this point. (See reference 25.) A somewhat different 
approach to this problem (reference 26) has led to a solution 
very like the one given here. 
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TABLE III.— SOLUTION OF EQUATION (122) 


>7 

■ 

/<») 


7 

- /M ’ 

/(si)*n 

0.0 

- 1.0000 

0.0000 

4.0 

-0.1307 

1.0124 

.2 

.9898 

.1990 

4.2 

-.1008 

.9893 

.4 

.9597 

.3940 

4.4 

-.0716 

.9720 

.0 - 

.9087 

.5808 

4.6 

-.0446 

.9604 

.8 

.8356 

.7552 

4.8 

-.0207 

.9539 

1.0 

.7339 

.9122 

5.0 

-.0008 

.9517 

l.a 

.0032 

1.0459 

5.2 

.0149 

.9531 

1.4 

.4597 

1.1622 

5.4 

■ .0285 

.9573 

l.fl 

.3188 

1.2300 

5.6 

.0340 

.9633 

1.8 

.1880 

1.2807 

5.8 

.0379 

.9705 ! 

2 0 

.0724 

1.3067 

6-0 

l .0388 

.9782 , 

2.2 

0245 

1.8115 

6.2 

.0373 

.9858 

2.4 . 

-.1008 

1,2900 

6.4 

.0338 

.9929 

2.8 

-.1662 ■ 

1.2733 

.. 0.6 

,0292 

.9992 

2.8 

-.1919 

1.2386 ■ 

6.3 

.0237 

I. 0045 

3.0 

-.2098 

1.1983 

7.0 

.0180 

1.0087 

3.2 

-.2126 

1.1561 

7.2 

.0124 

1.0117 ■ 

3.4 

-.2031 

1.1145 

7.4 

.0072 

1.0137 

3.0 

1843 

1.0758 

7.8 

.0026 

1.0146 

* 3.8 

-.1693 

1.0414 

7.8 

-.0012 

1.0148 



Fiopre 54.— Variation of fW with n 


necessary, therefore, to operate further on the solution given 
for the loading in the steady-state problem to obtain the 
solution for the loading in the unsteady problem. But 


(it, 


dAp" 

=-f 

1 r&wn 

_ bs 

i bsjo 

_ arj, 


so that £f the notation 

; 7 ;- : hrj. 


[~dA<p~] = Vo ( &P\ 
\?o/. 


is adopted (where ( — ) represents the loading in the anal- 
\9'o/» 

ogous steady-state problem), then the expression for the 
unsteady loading can be given in terms of (j~) by the equa- 
tion 

(%);''■] (124) 


By the application of equation (124) to equations (120), 
the loading for the various regions in the yt plane of the un- 
steady wing can be found. For region 2 in figure 55 there 
results 


^-R-,[T arctan Vf=v+w +m - u "b.(r 


The unsteady- analog, sinking wing. — The first step in 
deriving the unsteady-flow results for the sinking wing from 
the steady solution is to replace x with f. In equation (122) 
this corresponds to replacing i\ with r where r is equal to 
f/s (equations (116)). The second step is to replace w/V 0 
with —a and to rederive the expression for loading coeffi- 
cient since in the time-varying problem it is expressed in a 
somewhat different manner than in the steady-state analog. 
In the unsteady case, as the triangular wing moves tlirough 
a lixed reference plane the local span intersecting this plane 
grows as a function of time and. equation (12), which repre- 
sents the partial derivative with respect to time with x fixed, 
must be expanded to the form 

a p 2 r£>A(p“i _ 2 /r&A?~] , pA*q &s\ 

To ~ tWoL bf Jx - TWoVL at J, L as iu) 

where |p^~J aQ d indicate derivatives taken at 

constant $ and f, respectively. Since s is equal to m(z+AI Q t), 
5s/2>f equals mM 0 , and there results 


which becomes 


- 8 “ ["mMo-*/- — yr+m tan -J 8 — tW 
q a vMoL V »-\y\ V t — 

The loa'ding coefficient can be similarly derived in the other 
regions so that finally, for the regions shown in figure 55, 



Fiopre 55.— Regions used to express loading on section of unsteady, subsonltvodged, trl- 

angular wing. 


Ap_ 


7o 


(i23) 


In the steady-state problems an analog to the term 
involving [M is missing, and the loading coefficient is 

given entirely by an operation equivalent to • If is 


Region 1 


Region 2 


A p 4a 

q 0 ~M o 


(125a) 


A £ 

?o 



t—8+\y\ 

\y\ 



(125b) 
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Region 3 


<lo 


,re V(=rfe +a .' c ta “ Vr=v-y~ f) 


t—s—y 

*+y 


< 125c> 


For the interval f>2s equation (121) must be considered. 
By means of equation (124), the expression for the loading 
coefficient can be written 


ir 0 [ 4a/(r) V i_ (f ) 1 4a/ (« ) V i- (?) dti _ 


which becomes 


A/) 4a 
q <] J-/"o 


V 1_ vy 



<1261 

where /(r) is the solution to equation (122). Notice that for 
large r (when the loading has reached its steady state), /(r) 

is zero and I /(ri)dri is unity. (See table III.) Hence, the 
Jo 

loading is given by the equation 


so that — w(Vq in the steady-state equations (120) and (121) 
becomes SsjmV a . Since the loading coefficient is still given 
by equation (124), there results for the conversion of equa- 
tion (121) the expression 


Ap_ 4 ams 

y« 2 — y 2 

which is the steady-state value 'for a slender triangular wing 
(equation (112) when £’=1). 

It is now possible to derive the chord loading factor P 0 
as defined by equation (92a), thus 



Placing equations (125) and (126) in this expression, it is 
found that for 0<r<2 


and this can be reduced to the form 
A p_ 4 8 


(l-rni/ 0 r)^/(r)+ 


?0 mM a V a 


(128) 


As in the discussion of equations (127), it can be seen that 
equation (128) becomes for the steady state (t large) 

A p 4$ / 2a J — y* 


A p 4 6 

2o I T o \-\S J ~2/ / 


P 0 =2(2-r)+4mA/ o r (127a) 

and for r>2 

P,)=7r(l — m A/ 0 r)/(r)-i-2mAjfoir ( /(ri)dri (127b) 

Jo 

Since the values of P„ given by equations (127a) and (127b) 
were derived using different methods, their magnitudes at 
r=2 are not equal. The final curve for P 0 must be con- 
structed by fairing the solution for r<2 into that for r> 2. 
Figure 56 shows these results together with the final curve 
chosen (solid line). 

The unsteady-analog, pitching wing. — When the wing is 
pitching at a steady rate about its apex, the equation for 
the vertical induced velocity on the plan form is 

w K =— {x+A£j)6 


and this can be shown to agree' with the steady-state slender- 
wing results given in reference 27. 

It is now possible to derive the chord loading factor P t 
as defined by. equation (92b) 

<p = mM a V 0 f A p 

■ Vx 2s*0 J-,2 o dy 

Using equation (128), one finds for r>2 

P 1 =t{1 — mA/or)/(r)4-37n_l/ 0 rJ' /(ri)cfri (129a) 

and a similar analysis based on equation (120) yields for 
0<r<2 

Pi~^2-r+4mA/ 0 r-|mM 0 r J ) (129b) 
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As in the case for P 0 , the two equations for Pi do not join at 
r= 2 and the final curve must be constructed by fairing the 
solution for r<2 into that for r> 2 (see fig. 57). 



T 


FifiUKE 57.— Variation of Pi with r. 


Discussion. — It is now possible to assess the accuracy of 
the solution for very slender wings in the interval 0< r<2 
by comparing the values of loading given by slender wing 
theory (equations (125)) with the exact loadings given by 
equations (111), (113), and (114), and also by comparing 
the slender wing value of the chord loading factor P 0 (equa- 
tion (127 a)) with the exact value given by equation (117). 
It is apparent that the approximate solution differs from the 
exact only by a stretching factor in the f direction. Hence, 
if r is replaced by r//3, and m (note m is proportional to 
yjt) by mjS„ where /S a is given by equation (110), then equa- 
tions (125a), (b), and (c) are identical with equations (113), 
(111), and (114), respectively, and, of course, equation 
(127a) corresponds to equation (117). 

Tina rather remarkable, result can be enlarged upon from 
another viewpoint. Suppose .that in the steady-state analog 
problem the wing had been flying at some Mach number 
other than -y/2, say M,. The solution to such a new problem 
could be obtained from the old one merely by applying the 
Prandtl-Glauert correction, that is, by stretching all dis- 
tances in the x direction by the factor l/jS« where /V=|l — 
M, A \. Such a procedure would convert, for example, equa- 
tion (127a) to the form 

Pod«=2(2 — r/^)+4m ftAf, r/j8. 

Finally, if P 0 is adjusted so that P 0 =4 at r=0, there results 
Po = 2 (2 — r/5«) T 4mA/or 

which is exactly the answer given by equation (117). It is 
possible to simplify the statement of the.. procedure by simply 
remarking: The exact results for A pjq 0 or P 0 in the interval 
0<r//3,<2 can be obtained from the approximate results 
for a very slender wing by making an effective Mach number 


correction to the right-hand side of equation (125) or (127a), 
respectively. 

It is interesting to pursue this concept even further. 
Consider a spamvise section of a triangular wing as time 
increases from the starting impulse. The primary wave 
fronts emanating from either side pass across the section, 
forming the Mach lines in the steady-state rectangular-wing 
analogy. For very slender wings these lines make a 45° 
angle with the trace of the side edge and are used to divide 
the plan form into regions as in figure 58. . Now find tbc 
actual position of these primary wave fronts as they form a 
trace on the section in the yr plane. A Straightforward 
calculation shows that these lines actually make an anglo 
equal to arc tan l//9« with the trace of the side edges. Hence 
the effective Mach number which is used to correct the 
slender-wing results in the interval 0<r<2/S a is that which 
makes the Mach linos of the steady-state analogy coincide 
with the actual trace of the primary wave fronts. 


r— Actual trace of 



1NDICIAL LIFT AND PITCHING MOMENT ON VERY SLENDER TRIANGULAR 

WINGS 


The lift coefficient for the sinking whig is given in the 
notation introduced in equation (92a) by the equation 


f" It Podx 

mc 0 ,/o A* o 


where P 0 has been determined in the last section ns a function 
of r=f/s and c 0 is the root chord. Consider the situation 
at a certain fixed time and let the x coordinate in the above, 
formula be fixed in the wing. Then set 



8 

mc 0 


( 130 ) 
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and as before 

T 0 =——=mM^T 

Co 


(131) 


where r„ is the number of wing-chord lengths traveled. In 
this way it is possible to obtain the expression 

(132) 

and similarly 

(133) 


The equation for the lift and pitching-moment responses 
(.where again the pitching moment is taken about the apex) 
on a pitching wing are 


and 



The values of P n and P t were taken from curves similar to 
figures 56 and 57 in the last section (using the faired curves 
in the vicinity of r=2) and the results for the indicial lift 
and pitching moment in terms of r a , the number of chord 
lengths traveled, are shown in Part IV. 

PART IV— RESULTS AND DISCUSSION 

TWO-DIMENSIONAL RESULTS 

The methods presented in the previous sections have been 
used to calculate the indicial lift and moment curves for tw r o- 
dimensional wings flying at Mach numbers equal to 0, 
0.5, 0.8, 1.0, 1.2, and 2.0. These results are presented in 
figure 59. 

Figures 59 (a) and (b) show the variation of the indicial 
lift and moment on a sinking wing. The initial value of the 
lift is given by the expression ia/M 0 . When the free- 
stream Mach number goes to zero this expression still is 
valid, the initial value being a pulse of force that occurs at 
f=0. . The final values of the lift are simply the two-dimen- 
sional, steady-state results given by the PrandtI-Glauert 
rule. Figures 59 (c) and (d) show the lift and moment 
variation on a wing pitching about its leading edge. These 
functions were not computed for d/ o =0.5. The results are, 
of course, subject to the restrictions of linearized, com- 
pressible-flow theory and, for example, the calculated 
responses given in figure 59 for sonic speeds must be con- 
sidered as being outside the realm of validity within a few 
chord lengths of travel. In application to high-frequency 
oscillations, however, the initial portions of the indicial 
curves dominate the response characteristics of the airfoil 
and calculations near il/ 0 equal to one need not be invalid. 


More accurate values of the indicial curves are given in 
tables I and II for the subsonic Mach numbers 0.5 and 0.8. 
These tables can be extended to larger values of r 0 by means 
of equations (71) and (72). 




(a) LHton a sinking wfng. 

(b) Pitching moment (about leading edge) on a sinking wing. 

Figcm 59. — Variation, of two-dimensional indicial lift and pitching-moment coefficients with 
chord lengths traveled for* evaral Mach numbers. 
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(c) Lift on wing pitching about Its leading edge. 

(d) Pitching moment (about loading edge) on a wing pitching about Its leading edge. 
FimjRE 59.— Concluded. 


TRIANGULAR WINGS WITH SUPERSONIC EDGES 

Curves of the indicial lift and moment on sinking and 
pitching triangular wings with all supersonic edges are shown 
in figures 60 and 61 for Mach numbers 1.2 and 2.0, respec- 
tively. For the purpose of comparison, the curves for a 
two-dimensional wing flying at the same Mach number 
havo been included in each figure, as well as the curves for 
the same triangular wing in reversed flow. 

Several conclusions can be drawn from these results. 
First, notice that, the total indicial lift on the triangular 
sinking wing is the same at every instant as that on' the same 
wing in reversed flow (both wings, of course, having started 
with the same velocity at the same time), and that the value 
of this lift is the same 'as the total indicial lift on the two- 


dimensional wing only at the beginning of the motion and 
again 'when the steady state has been attained. 

Secpnd, notice that, since all the characteristics for the 
triangular wings are independent of the nnglo of sweep, they 
are vft Kd for any unyawed triangular wing as long as the 
edges are supersonic. 

Third, it is apparent that the transition of the total indicial 
lift frpm_its initial to its final value is less abrupt for the 
triangular than for the two-dimensional wing. The initial 
. and final values of C L depend on l/M 0 and 1//3, respectively, 
1 so that as the Mach number is increased the variation dies 
out altogether since ft and M 0 approach one another. The 
same remark applies to all the other coefficients. 

' " SLENDER TRIANGULAR WINGS 

The results for the lift and moment on a slender triangular 
wing that is sinking or pitching about its apex are shown in 
figure, 62. The analysis by which the results were obtained 
is valid when both rn, the tangent of the semiapex angle, and 
mM 0 are small; the results are given for 1/8.' The 

curved &re all qualitatively alike, in each case the response 



(a) Lilt on »t nk[n g wings. 

Fioube 60.— Indicia] aerodynamic characteristics of triangular wings with supersonic edges; 

IWi-1.2. 
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(c) Lift on wings pitching about their leading edge or apex. 


Fioore 60.— Concluded. 


(d) Pitching moment (about leading edge or apex) on wings pitching 
about their leading edge or apex. 
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0 .4 .8 12 l.e 80 2.4 

Chord lenqfhs traveled, V a t'jc 

(c) Lift on wtn*s pitching about their leading edge or apex. 

(d) Pitching moment (about leading edge or apex) on wings pitching about their leading 

edge or apei. 

FiGtns 81.— Concluded. 
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.5 12 1.6 2.0 24 

Chard lengths traveled, V„tfc 


Chord lengths traveled, V a t'jc 

(Si Lift, 

fbi Pitching moment (about leading edge), A&»U. 

Ftotm* 63.— Indlcial lift rind moment on sinking rectangular wings flying at supersonic 

«p#wis. 

RECTANGULAR WINGS 

By means of the analysis and curves presented by J. W. 
Miles in reference 2, the indicial results for a two-dimensional, 
sinking wing can be compared with those for a rectangular, 
sinking wing flying at supersonic, speeds. This comparison 
is given in figure 63 for Mach numbers equal to 12 and 2.0 
and aspect ratios equal to 2, 4, 6, and infinity. It can be 
seen that the modification in the shape of the indicial curves 



Chord lengths traveled, V.t'jc 


(c) Lift, Af»-2A 

(dj Pitching moment (shout leading edge), A/i—2.0. 

Fioubi 83.— Concluded. 

brought about by the decrease in aspect ratio is similar to 
that caused by a decrease in Mach number from supersonic 
to subsonic magnitudes. 


Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., October IS, 1950. 
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APPENDIX A 


DETERMINATION OF SUBSONIC, INDICIAL, 
SECTION LIFT AND PITCHING -MOMENT CURVES 


THE LOAD DISTRIBUTION 

The following results for the indicial load distribution on 
sinking or pitching wings can be obtained in two ways. 
One of these methods will be outlined in the subsequent 
paragraphs. The other is outlined in references 28 and 29 
and is referred to as the lift-cancellation technique. The 
latter method has been used to check the load distributions 
originally obtained by the former so that an independent 
check of these results has been carried out. 

It was shown in Part II, Two-Dimensional Boundary- 
Value Problems, that the lifting-surface analog to the solu- 
tion for load distribution over an unsteady, two-dimensional 
wing traveling at a constant subsonic speed involved the 
calculation of load distribution over a swept-forward wing 
tip with subsonic edges. Figure A1 indicates the geometry 
associated with the boundary conditions. Solutions are 
given only for the loading in the five regions shown. As has 
been discussed in the text, for il/ 0 =0.8 this was considered 
adequate to define the behavior of the indicial responses in 
the early stage of motion. For M a — 0.5 the loading in regions 
farther down the wing had to be calculated; and in all, for 
this Mach number, the loading was analyzed in the 11 
regions shown in figure 25. However, the latter analysis 
was carried out by a numerical application of the lift- 
cancellation technique and none of the details will be pre- 
sented here. 

In the notation of the unsteady problem the expression 
for the velocity potential can be written 


1 C C w^dtjdxi 

J S-uy-ix-xtf 


(Al) 


where r is the area on the wing plan form included in the 
Mach forecone from the point ( t,x ). Equation (Al) is 
applicable only for cases in which »„ is known at all points 
within the forccone, as is the case when the edge of the wing 
within the forecone is everywhere supersonic (i. e., region 1 
in fig. Al). However, Evvard (reference 30) has extended 
the solution provided by equation (Al) to include cases such 



Finmtjc Al.— Regions used In analysis of unsteady, two-dimensional wing flying at subsonic 

speed. 


as shown hi figures A2 and A3 hi which the forecone inter- 
sects a subsonic edge and includes a region of unknown 
upwash. As was pointed out in reference 30, equation (Al) 
applies in these instances if the area of integration t is limited 
to the shaded regions shown in the figures. It is apparent, 
therefore, that the potential (and thus the loading) over a 
sinking or pitching wing can readily be determined for regions 
1, 2, and 4, hi figure Al. 

Points in regions 3 and. 5 have forecones which intersect 
two subsonic edges, and the method just discussed can no 
longer be directly applied. In reference 30, however, a 
method was given of evaluating the upwash in the region 
between the Mach cone from the apex and the leading edge 
(region 6 in fig. A4). Thus, the plan form has become, 



Fioube AS.— Axes ol Integration when Mach foreconc Intersects subsonic trailing edge. 



Fiocre A3.— Area of integration when Mach forecone Intersects subsonic leading edge. 



Fracai A4.— Region 11, area between wing leading edge and trace of foremost Mach gone. 
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effectively, one such as shown in figure AS in which only one 
edge is subsonic. This reduces the problem of finding the 
potential in these regions to the same problem as was involved 
in region 4. The analysis used in finding the loading over the 
various regions will now be considered. 


wing and pitching wing, respectively. The method developed 
in reference 30 was used to obtain the value of w u over the 
area between the lines a=0 and r=«(l+4/ 0 )/(l— AQ (region 
6 in fig. A4). The results for the sinking and pitching wing 
are respectively, 



Fuwbe AS —Region of Integration for mired sonic and supersonic leading edge and subsonic 

trailing edge. 



.r 

Flora* Ad— Regions of figure A I in the rt plane. 


First, introduce a new coordinate system in which the lines 
x=—t and x—t are taken as the r and a axes, respectively. 
(See fig. A6). This amounts to a rotation of the original 
system of axes through an angle of 45°. The transfor- 
mations relating the r,s to the x,y system are 

r=\(t-x) t~(r+s) 

s I v l 

a=4=(f+ar) x==-^(a— r) 

In the new coordinate system equation (Al) is written 

1 f C wujru&ddnds, 

t^2Jt J V(r— rO(a— «0 


2 tv r / 2 s 

W ’~ t L V (r—s)—M 


MoCr+a) 


arc tan 


/ 2a 

\(r— a)— Mo(r+a)_ 


w, 




M^r+s) 


[0-„)-M„(r+,)] arc t« - V / (r _, ) _V f + ,j 


(A3) 


In the r,a coordinate system, if the subsonic trailing edge is 
not included in the forecone from the point at which the loading 
is being determined, the expression for the loading coefficient 
can be written 

A? -2 /d_ d\ f f igufri.aQcMa, f . . 

(Zo irT r oA/oV^ r da/J r J y(i — ri)(a— «i) 

However, if the trailing edge is included in the forecone, 
and if the Kutta condition is to apply along such an edge, it 
can be shown that the equation for the loading coefficient 
assumes the form 


A V 2 r dr (b . d \ f- iP.Crt, aQffot . 

2o ir T’VfoJry'r — r t v^ r i sjs — 8i 

SIN KIN G WITS G 

The preceding method can be applied to the sinking wing 
to obtain the following integral relationship for the loading 

over the various regions. The subscripts on (~^ indicate 

the region for which the particular equation applies. 

1_/Ap\ = J^/jd d\ f r dr ' f* dSl 

a \ 2o / 1 ir-TIp da/ — ri J- r -\8— 

i / Ap\ 2 /a , s\ r» d Sl r dr , 

ff\ 2 o/i wMo\dr da/ J rq— Aft) ya — a, J -*t V — r t 

(l-t-Aft) 

1 /A £\ 2_ ( r dr j /A + A) f* Jh ■ 

a V 2o/» ffil/o ( J -yjp — \dr i da/ J nq-Aft) ya— a. 


(l+Aft) 


nq-Aft) 
1 0+355 


-ds, 


Jo y'r — r t i da/ J 

I 2 a t / 2e t 

yjri-l i)-^oOvtgi) arC an V (ri— ai)— MoC/yhaQ 


\ ! s—s i 


+ 


ro dr t / d d\ C* <fa t ) 
)r d v’r— r x l ^8/ J-ri y 1 , 8—8i ) 


The vertical induced velocity e, over the wing plan form 
is given in equations (2) and (3) of the text for the sinking 
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« \<?o/ 


1 \ — 2 f r & r >. ( f?__L <■> \ f* 
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where 


/ 28 , 


/ 2Sl i 

V (ri-Si)-M 0 (ri+«i) 

ul u tmi'i 

/ (ri— * 0 — M 0 (r, 4 -Sij } 

_ 

V® ®i 
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1 — M 0 


Most of these integrals can he readily evaluated to give 

(A6a) 


1 /Ap\ = J_ 

«\2o/i M 0 


1 = 8 (-M* - f 
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where 
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4= arc sin V-To+ Af 0 4 
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PITCHING WING 


arc sin r(r+®Oa- Mo)[(8--8 1 )(l+M 0 )-cV2l+28 1 [8(l + Af 0 )- C ^/2-Kl-A/ 0 ); 


A similar analysis of the pitching wing yields the following 
results for the loading coefficient: 

(A7a) 

q\q o/a 3 tA1qL(14-A/o) VMo4+:Eo 

3(Mo4+*o> arc tan l M <t°+ x ° (A7b) 
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LIFT AND PITCHING-MOMENT COEFFICIENTS 

The lift and pitching-moment coefficients may be obtained 
by suitable integrations of equations (A6) and (A7) and 
are given in the time intervals indicated in figure 34 by the 
following expressions: 

Sinking wing 

0<4< 1 


1 +44 


[l-f 0 (l-A/ a )I 
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Pitching wing 
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APPENDIX B 


DETERMINATION OF SUPERSONIC, INDICIAL, 
SECTION LIFT AND PITCHING -MOMENT CURVES 

THE LOAD DISTRIBUTION 

In the case of the unsteady supersonic wing the expression 
for the velocity potential may be readily obtained by placing 
the values of w v given by equations (2) and (3) of the text 
in the equation 

-ff w - ^ Xl (Bl) 

where r is the area on the plan form included in the Mach 
forecone. The loading may then be calculated from the 
relationship given in equation (12). 

SINKING WING 

The load distribution over the regions A, B, and c shown 
in figure 34 are given by the following expressions: 
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PITCHING WING 

For the case of the pitching wing the values of ^ (ji~) 
in regions A, B, and C are, respectively, 
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LIFT AND PITCHING-MOMENT COEFFICIENTS 

The lift and pitching-moment coefficients may be obtained 
by suitable integrations of equations (B2) and (B3) and 
are given in the time intervals indicated in figure 34 by the 
following expressions: 
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